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ABSTRACT 


This thesis is concerned with the use of computer 
generated binary masks, invented by A. W. Lohman, as matched 


filters. 


A theoretical analysis of the transmission 
properties of a filter previously invented by A. B. Vander 
Lught is presented and it is shown that the Lohman filter 
can be regarded as a binary simulation of a Vander Lught 


filter. 


Using fingerprints as test objects, the 
computational requirements for the evaluation of parameters 
describing a Lohman filter are examined and the production 


of a filter is explained in detail. 


Several Lohman filters designed to detect 
fingerprints are made and experimentally tested. The 
experimental results demcnstrate the ability of a Lohman 
mask to identify fingerprints and show that the accuracy 
achievable with the experimental system used in the present 


work permits 49 prints tc be processed in parallel. 
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CHAPTER 1 


1.1 Spatial Filtering 

In recent years the interest in optical data 
processing, or optical ccmputing as it has come to be 
called, has increased rapidly. Amongst the most successful 
branches of optical computing is the field of spatial 


filtering [1]. 


The meaning of spatial filtering will be explained 
later in some detail but at this moment it is sufficient to 
say that a spatial filter is a two-dimensional transparency 
or mask which in the most general case is required to 
control, in a predefined manner, the spatial distribution of 
both the amplitude and phase of a light wave passing through 
oe os The attenuation of amplitude can be easily 
accomplished by the use cf, for example, a darkened 
photographic plate. The control of the phase of a light 
wave was, however, a very difficult task until holographic 


techniques were developed. 


The spatial filter considered in this thesis is a 
type of a matched filter. Historically, the matched filter 
was developed in electronic communication as the filter 
which maximizes the signal to noise ratio when detecting a 
Signal submerged in white noise [2]. The usefulness of an 


optical matched filter is due to other reasons which follow. 
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The Fourier transform plays an important role in 
the theory of matched filtering. While the Fourier 
transformaticn of electric signals is a rather complicated 
operation, requiring extensive hardware, a single lens can 
perform the Fourier transformation on any optical signal, 
that is the two-dimensional distribution of magnitude and 
phase of the optical disturbance. Ttees ths ability ofea 
lens to readily perform a Fourier transformation which makes 
an optical matched filter an attractive tool regardless of 
the statistical distributions of the noise in which the 


Signals to be detected are buried. 


Prior to holographic techniques, spatial filters, 
including matched filters, were generally produced in two 
separate stages. The attenuation of amplitude was 
accomplished by the use of a darkened photographic plate and 
the control of phase was generally achieved by the use of a 
transparent film with appropriate variations in thickness. 
With these techniques it was only possible to implement 


rather simple transmission functions. 


In 1963, A.B. Vander Lught invented a holcgraphic 
method of producing masks consisting only of an absorption 
pattern [3]. These masks therefore only affect the 
amplitude of light passing through them but nevertheless 
Vander Lught demonstrated that they could be used as matched 
filters despite the fact that they do not affect phase. 


This technique is therefore a method of reducing the 


aso anal ekpaie » qoxbwored wrivinoae vation 40 Bfisa6go 
loaphe Lebiitqe pte co nvtyseso%ne7? TetTiot sat azotzeg 
bac shvsiapaw Xo noktudeatelt Ledorenseto-ows sat ot taue 

5 20 YSttida cide at 9T ~epedaegarh ined tao d+ Yo sastg 
seien i9hde aokiaezeteweds tsisifet 5 mtotteq ylibso2 oF eoat 
160 seelboaped /teos evisueitin né yeoilt baiiotan Lsoitqo me 


ait? worse az eetat sit 30 enottudriteih Lsotteionte eis 
a 


7 
f 


f belied st6 fetooteh oi of aleapea 
j 

,etertea. Lebteye ,ase— tadnet baemnne tod of 20629 
ow! at botiubeth Vileaetsp aie ,2magini fadorse eniavisuk 
soy sbUeitjae 2 moiYausesss si? 8 .espeta sreTeyee 
bus stelq Siduarpodosy Betedtah 6 2 seh edd YA bade ciquesms 
6 20 eeu ait yd beveltips yilavedep sev Ssesty Yo Iottios ede 
Jeaeada> il? ni enoigataey esstiqowige a2iv siti sve te_anE7? 
dee ge: oy sidgeacg yluo abe 4% coupindses cases ieee 
2002 70a0? dodeadsensit edqgats tew?rer 


skdgewpo log @ betaevai tiqwd ashaey of. ,F Aer at 
HOLTELOOaN ae fo ¥ino gaivelenos 2deee pabothocy 46 bottes 
‘ed? JosiMbs (iyo Stotexsd? adaow dest? [©] azesenq 
sasibisrbena tue weld dpyordt pT isesy riplt Io sbuaLEgee 
Pototee @b Neeu od Diwoo yald ssd9 Godeadeqo met Sanwa TobaDY | 
«seedy teks 200 ob ysis Frdd tIet 907 EItGees sxoStht 
ode patsobex to bodsen 8 az0%0400) ot cuphedse® alAt 


generally complex transmission function of a matched filter 


to a real transmission function. 


The Vander Lught method was the first convenient 
method of producing optical matched filters and requires 
that the signal to be detected (the signal to which the 
filter is matched) is available in the form of a 
photographic transparency. This restricts the use of the 
Vander Lught technigue tc real Signals, that is signals 


which consist of amplitude distributions alone. 


The Vander Lught filters have been applied to 
various signal detection problems including the recognition 
of alphanumeric characters [4] and fingerprint 
identification [5,6]. The experiments with fingerprint 
identification have shown particularly interesting results. 
These experiments show that despite the complexity of 
fingerprints they are detected with a high degree of 
accuracy by the Vander Lught filter. Even if only a 
fraction of a fingerprint is available to construct the 
filter, recognition is pcessible [5]. Also, the filter 
indicates similarity when the print to be recognized is 
superimposed by a second print [6]. In the last case the 
print to be detected was so badly obscured by the 
superimposed print that it could not be identified by visual 


observation. 


These properties of the Vander Lught filter are 


very important if they are to be used in criminal 
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investigations, since the prints collected at the scene of a 
crime are often incomplete or obscured by superimposed 


prints. 


Systems have been developed which use the Vander 
Lught filter to identify a particular print from a file of a 
large number of prints [7]. In one experiment where such a 
system was used to identify 64 prints from a file of 10,000, 
in 65% of the cases the correct print only was identified 
and in 90% of the cases the system reduced the number of 
prints that had to be manually searched to less than 0.24% 


of the complete file. 


In 1966, A.W. Lohman reported a method of 
producing optical filters (masks) consisting of only two 
different transmission levels [8]. These masks were 
originally intended for use as holograms but the inventors 


Showed that they could also be used as matched filters. 


This method of making matched filters uses a 
digital computer to evaluate parameters describing the 
required filter. The filter is then plotted by a computer- 
guided plotter in a greatly enlarged scale and finally this 


plot is photcgraphically reduced to the desired size. 


Some of the applications which the Lohman mask has 
been used for are holography [9], code translation [10] and 
various filtering operations [11] including edge 


enhancement. 
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The main advantage of a Vander Lught filter is 
that when the signal to he detected is available on a 
photographic transparency the filter is very easily 
produced. In the experiments mentioned previously the 
Signals to be detected (letters and fingerprints) were 
available and the production of a Vander Lught filter was 
straightforward. For this reason rather limited attention 
has been paid to the Lohman masks as matched filters and no 


comprehensive study of their properties has been published. 


In many respects, however, the Lohman filter has 
advantages over the Vander Lught filter. In order to 
produce a Lohman filter it is only necessary that the signal 
is known in mathematical terms. This is particularly 
important if the signal to be matched is not a real function 
but a complex one (i.e. the signal is both an amplitude and 
phase function). In this case a Lohman filter can easily 
be produced whereas the difficulty in obtaining the complex 


signal in the form of a transparency is considerable. 


The situation described above, where it is desired 
to detect complex signals, frequently arises when optical 


components are to be tested [12]. 


If the signal to be detected is known to be 
submerged in noise, other than white noise, a filter which 


is optimum in the sense that it maximizes the signal to 
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noise ratio, can easily be made by the Lohman technigue. 
All that is required is that the non-white noise 
distribution is allowed for when the filter is calculated. 
Such an optimum filter can also be produced by a 
modification of the Vander Lught technique {3] but the 
procedure is much more cumbersome than when the Lohman 


method is used. 


It may be required to produce a filter matched to 
a Signal which is not actually known but is to be found 
through experimentation or some sort of a learning process. 
In this case a Lohman filter could be synthesized directly 
whereas to synthesize a Vander Lught filter a signal would 


first have tc be selected and realized. 


The advantages of the Lohman filter described 
above only apply to its use as a matched filter, but it 
should be reiterated that the Lohman mask also has many 


other applications as indicated in the previous section. 


There are therefore many instances where a Lohman 
mask is more Suitable than a Vander Lught filter. For this 
reason it is desirable te investigate in some detail the 
performance cf the Lohman filter and compare it to the 


Vander Lught filter. 


1.3 Scope of the Thesis 
The scope of this thesis is to study the 


relationship between the two filters previously mentioned. 
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In Chapter 2 some theoretical background is given so that 
the content cf the remaining chapters may be better 
understood. In Chapter 3 a detailed description is given 
of both the Lohman and Vander Lught filters and it is 
explained how a mask consisting only of an absorption 


pattern can act as a matched filter. 


Also in Chapter 3 the relationship between the 
Vander Lught and Lohman filters is examined, and it will be 
shown that the transmission functions of the two filters are 
so Similar that the Lohman filter can be regarded as a 


binary approximation to the Vander Lught filter. 


When a Lohman filter is being produced, the signal 
must be sampled in order to enable the computing of the 
filter parameters. This leads to the question of how many 
Samples are necessary to adequately represent the signal and 
what the computing costs involved are. This problem is 
considered in Chapter 4. As mentioned previously the 
results of experiments with fingerprints have shown 
promising results so all the experiments performed are done 


using fingerprints as test objects. 


In Chapter 5 the various steps in the production 
of a Lohman filter, including the computation of the filter 
parameters and plotting cf the filter, are explained in some 
detail and in Chapter 6 the results of fingerprint detection 
by Lohman filters are presented. The results of these 


experiments are compared with the published results of 
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similar experiments carried out with Vander Lught filters. 


Finally, in Chapter 6 both theoretical and 


experimental results are summarized and discussed. 
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CHAPTER 2 


2.1 The Fourier Transfors 


As previously mentioned the Fourier transform is 


an important concept in the theory of spatial filtering. 


If g(X,y) is an optical signal in a flane with 
Spatial coordinates (x,y) the Fourier transform of g(x,y) is 


defined by 


F{g (xey)}=f[.9 (x ¥) Exp[-j2(uxtvy) Jaxdy (2.1) 


where F{ } denotes "Fourier transform of". A Fourier 
transform is also referred to as a Fourier spectrun, 
frequency spectrum or spectral density. The coordinates u 
and v in the Fourier transform domain are referred to as 


Spatial frequency coordinates. 


The Fourier transform has an inverse which is 
defined, for any complex function G of two independant 


variables u and v, by 
F-1{G (u,v)}=/f_G (u,v) Exp[ j2m7(ux+vy) Jdudv (2.2) 
where F-1{ } denotes "inverse Fourier transform of". 


There exist several theorems expressing certain 
properties of Fourier transfcrms. Some of the most 


important of these theorems are: 
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Te Linearity Theoren. 
If g and h are two complex functions of two independent 
variables x and y and a and b are complex constants 


then 


F {axg (x,y) +bxh (x,y) }=aF {g (x,y) } +bF {h (x,y) } (2.3) 


Fle Convolution Theorem. 
If F{g(x,y)}=G(u,v) and 
F {h (x,y) }=H (u,v) then 


F{f {9 (Ped) h (x-p, Y-9) dpdq} =G (u,v) H (u,v) (2.4) 


The quantity inside the curly brackets on the left hand 
Side cf the equation is, by definition, the convolution 
of the functions g and h. 

36 Fourier integral theoren. 
If g(x,y) is a continuous function then 


F{F~-1{g (x,y) }} =F-*{F{g (x,y) }} =9 (X-Y) (2.5) 


A statement of these theorems and others together 


with proofs is given in reference [13]. 


2.2 Linear Systems 


In optical data processing one is frequently 
concerned with the relationship between a two-dimensional 
complex amplitude distribution (sometimes referred to simply 
as "light distribution") in the input plane of an optical 


system and that of the outfut plane of the systen. 


Any optical system can be represented as an 


operator P which acts on the input to produce the output. 
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Thus if 1(X,«Y,) is the light distribution in the input 
plane and o(x,,y,) is the light distribution in the output 


plane, the effect of the system can be expressed as 
O(X5¢Y5)=P {i (x,+Y,)} (2.6) 


As an example i(Xj0Y,) could be the light distributicn in 
the object plane cf an inaging system and O(X5,y,) the light 


distribution in the image plane. 


A system is said tc be linear if the following 
expression is true for all complex constants a and EF and all 


inputs A(x 4.74) and J(%,°Y,): 
Pfaxi (Xj e¥,)+b*5 (XY ,)} 


“Par {a (xyey i} +bP (d(x pe¥ 23 (2.7) 


The linearity cf a system implies that any input 
can be decomposed into as many parts as desired provided 


only that the parts add up to the actual input. 


The Dirac delta function is a convenient way of 


decomposing an input. The Dirac delta functicn may be 
defined by 
6(x,y)=Lim N@Exp[-N271(x2@+ty2) ] (2.8) 
Noo 


Making use of the Dirac delta function the input i(x,-Y,) to 


an optical system can be rewritten as 
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i(xy-y p=! l.4 (pea) § (x y-p- 1-4) apag (2.9) 


The property of the delta function which makes it so useful 
for the purpcse of decomfosition is called the sifting 


property and is discussed in detail in reference [13]. 


Combining equations (2.6) and (2.9) the output of 


a linear optical system can be written as 


O(XoeYo)=ff i (Ped) P {6 (x y-P-¥ 4-4) }dpdg (2.10) 


In the above equation the linearity property has been used 
to bring the operator P under the integral. It should be 
pointed out that while the right hand side of equation 

(2.10) does not depend explicitly on x» and y», the result 
of the operation by P is a function of x» and y» as can he 


seen from equation (2.6). 


Digs quantity PAc(Xy—-p,yj-d)} is the output Of the 
system due to an impulse (delta function) at coordinates 
(p,q) in the input plane. This is called the impulse 


response of the system and will be denoted by h. Thus 
h(X oe Y 5¢Ped) =P {5 (X 1-P-¥ 1-9) } (2.11) 


A system is said to be space-invariant if the 
impulse response does not depend explicitly upon all the 
coordinateS Xo, Yo, p and q but only on the differences 


(Xo-p) and (yo-q)- In this case equation (2.10) becomes 
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o(x ,6¥,)=f [2 (Ped) h (x,-pey ,-4) dpag (2.12) 


Comparing equation (2.12) to equation (2.4), the right hand 
side of equation (2.12) is seen to be the convolution of i 
and h. From now on the convolution operation will be 


denoted by * and equation (2.12) can be rewritten as 
O(X,ry,)=i*h (2.13) 


Using the convolution theorem (2.4) equation (2.12) can be 


Fourier transformed: 
O(u,v)=H (u,v) TI (u,v) (2.14) 


The upper case letters denote the Fourier transform of the 


functions indicated by the ccrresponding lower case letters. 


Equation (2.12) is a space domain description of 
the system while (2.14) is a frequency domain descriftion. 
H(u,v) is the Fourier transfcrm of h(p,q), the system's 
impulse response, and is called the transfer function of the 


systen. 


The computational task of calculating a Fourier 
transform is a rather complicated one for all but the 
Simplest functions. However, a very Simple optical system 
consisting only of a lens and free space can perform the 
two-dimensional Fourier transformation on any light 


distribution at a very high speed [ 14}. (Essentially the 
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time taken tc perform the transformation is the time taken 
by light to travel from the input to the output of the 
optical system). This property of an optical lens makes it 
an extremely powerful tool for the purpose of optical 


Spatial filtering. 


Figure 1 shows an optical configuration which can 
be used for any spatial filtering task, including matched 


filtering. 


Plane P1 contains S(X,0¥,)s the input to the 
system, i.e. the pattern that is to be filtered. S(X4/Y}) 
is illuminated by a plane wave of amplitude a, so the 
complex light amplitude distribution b(x,+Y 3) of the wave 


emerging from P1 is given by 
b(X,/Y,) =axs(Xi+Y,) (2.15) 


The planes P1 and P2 are each separated from the 
plane of the (thin) lens by f, the focal length of the lens. 
Under these circumstances L1 performs a two-dimensional 
Fourier transformation on b(X,/Y,) {14] and the Fourier 
spectrum of b(x,+Y)) is displayed in plane P2. The light 


amplitude C(X,rY,) at plane P2 is given by 


¢(X5/Y,)=_a_S (u,v) (2.16) 
jrt 


where S(u,v) is the Fourier transform of S(x,,/Y,) [ 14}. 


The spatial frequency coordinates (u,v) in plane 
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P2 are related to the real coordinates (X50Y >) by 
u=x,/ (Af) and v=y,/ (Af) (Zt7) 


where A is the wavelenght of the coherent plane wave which 


illuminates S(X,°Y,)- 


Whenever the light distribution in one plane of an 
optical system has a Fourier transform relationship with 
that of ancther plane of the system it is convenient to 
consider one of the planes in terms of the frequency 
coordinates. This is frequently done throughout the 
remainder of this thesis. Under these circumstances it 
must always be born in mind that the two types of 
coordinates are related Ey equation (2.17). The unit of 
the real coordinates is that of length whereas the unit of 
the spatial frequency cocrdinates is 1/(length). 
Coordinate-axis in figures in the thesis are labelled with 


the coordinates most freguently used in the text. 


In spatial filtering the absolute amplitude of the 
Fourier transform in plane P2 is of little interest and thus 
the constant a/(jAf) in equation (2.16) is frequently 
dropped. Throughout the remainder of this thesis the 
constant will be neglected and thus C (Xo Y 5) is assumed to 


be given by 


c(x )=S (u,v) (2.18) 
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A filter is introduced into plane P2 of Figure 1, 
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thereby altering the frequency content of s(x 1,Y}1)- This 
modified spectrum is then Fourier transformed once more by 
lens L2 and the filtered version of S(Xj/Y ,) appears in 


plane P3. 


The lens L2 performs a Fourier transformation 
whereas an inverse Fourier transform is reguired. 
Comparing equations (2.1) and (2.2) it is seen, however, 
that an inverse Fourier transform is in fact a Fourier 
transform in reflected coordinates. (X is replaced by -x 
and y by -y). This is the reason for the reflecticn of 


coordinates (X33) in plane P3. 


It should be pointed out that the transfer 
function of this optical system is simply the transmission 


function of the filter placed in plane P2. 


A filtering system is said to be matched toa 
signal s(x,y) if it has a transfer function given by S*(u,v) 
[2]. S(u,v) is the Fourier transform of the signal s(x,y) 
and * denotes a complex conjugate. Since the transfer 
function of the system shown in Figure 1 is the transmission 
function of the filter in plane P2, the above defintion 
implies that a filter matched to s(x,y) must have 


transmission function S* (u,v). 


When a filter with transmission function S*(u, v) 


is present in plane P2 and s(x,y) is present at the input, 
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the light amplitude d(x5,y>5) behind plane P2 is given by 
d(X5¢Y5)=S (u,v) S* (u,v) (2.19) 


The right hand side of equation (2.19) is necessarily a real 
quantity. This implies that the light amplitude 

distribution behind plane P2 has a constant phase throughout 
the plane i-e. is a plane wave. A plane wave emerging from 


plane P2 is focused to a point in plane P3. 


If any transmission function other than s(x,y) is 
present at the input, equation (2.19) is no longer true and 
the wave emerging from plane P2 is not a plane wave, and 
will not be focused to a point in plane P3. Since a spot 
of light in plane P3 indicates that s(x,y) is present at the 
input or that s(x,y) “has been recognized" at the input, it 
is referred to as a recognition spot. The above also 


explains why the filter is called a matched filter. 


A matched filter can also be defined in terms of 
its impulse response. It is assumed that the filtering 
system is sapce~invariant and thus as can be seen from 
section 2.2 the impulse response of the system is given by 


the inverse Fourier transform of the transfer function. 


It is a well known property of Fourier transforms 


[15] that 


Fist (=x ,7y)}=5" (u,v) (2.20) 
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where S(u,v) is the Fourier transform of s(x,y). A matched 
filter is then required to have an impulse response given by 


S'(*x,~y)- 


Consider now the optical system shown in Figure 2. 
Lens L1 Fourier transforms the spherical wave diverging from 
the point source in plane P1, so plane P2 is illuminated by 
a plane wave. Lens L2 performs an inverse Fourier 


transform on the light emerging from plane P2. 


A matched filter can now be equivalently defined 
in terms of the optical system of Figure 2 as any mask which 
when introduced into plane P2 gives rise to s*(-x,-y) 
appearing in plane P3. Two comments should be made about 
the last definition of a matched filter. First, it ic 
strictly required that s'(-x,-y) and only s'(~x,-y) appears 
in plane P2 when a filter matched to s(x,y) is present in 
plane P2. However, if s'(-x,-y) appears somewhere in P3 it 
is always possible to shield off any light distribution 
which may appear elsewhere in the plane. Secondly, since 
the system of Figure 2 is assumed to be space-invariant, the 
point source in plane P1 may be displaced off the optical 
axis. The only effect of this is to tilt the plane wave 
incident on F2 and to displace the complete light 


distribution in plane P3. 
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2.6 Production of Matched Filters 

In order to realize an optical matched filter it 
is required to produce a transparency with transmission 
function given by S'"(u,v). This is in general a complex 


function, i.e. it is required to control both the amplitude 


and phase of light passing through the filter. 


While S*(u,v) is a complex function, S'+tS is 
always a real function. Technigues exist by which S*+S are 
recorded on photographic film in such a way that when the 
film is inserted into plane P2 of Figure 2, the responses in 
plane P3 due to the two terms S*’ and S are separated. Two 
of these techniques were outlined in Chapter 1, and will be 


examined in mere detail in the next chapter. 
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CHAPTER 3 


Vander Lught filter is a photographic record of 
the interference fringes due to a light amplitude S* (u,v) 


and a tilted plane wave. 


Figure 3 shows a schematic diagram of an optical 


System which could be used to record a Vander Lught filter. 


A mere practical configuration is the modified 
Mach-Zehnder interferometer criginally used by Vander Lught 
(3] and other configurations later shown by Vander Lught 


[16]. 


Referring to Figure 3, plane P1 contains a 
transparency with a transmission function S'(-x,-y). This 
is the impulse response which the filter is required to have 
and is generally a real function. The lens performs a 
Fourier transformation on s'(-x,~-y) and displays S* (u,v) in 


plane P2, the back focal plane of L. 


The part of the plane wave which passes above the 
mask in P1 is deflected by means of the prism to strike the 
plane P2 at an angle 6. The complex amplitude resulting in 
plane P2 from this tilted plane wave (sometimes called a 


reference wave) is given by 
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U(X) ¢Y,) =aExp(-j27wx,) (S571) 
where a is the amplitude of the plane wave and w is given by 
bess) 

w= (Sind) / (3-2) 


and i} is the wavelenght cf light used. 


To produce a Vander Lught filter, a photographic 
plate is inserted into plane P2 and exposed to the light 


amplitudes previously described. 


It is assumed that after development, the plate 
will have transmittance proportional to the intensity of 
light which expcsed the flate. The transmission function 


of the finished filter can then be written as 
T (u,v) = abxp (~j2mwx 1) +S" (u,v) |? (3.3) 


In the above expression a constant of proportionality is 
neglected. Expanding equation (3.3) and replacing x, by 


uvAf yields 


T(u,Vv)=a244S* (u,v) | 2+aS* (u,v) Exp(j27w\fu) 


+aS(u,v) Exp(-j27w\fu) (3.4) 


If the complex function S* (u,v) is separated into amplitude 


and phase by 


S* (u,v) =A(u,v) Exp(-j@ (u,v) ] (3.5) 
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the transmission function can be rewritten as 
T (u,v) =a2tA2 (u,v) +2aA (u,v) Cos[ 27wAfu-G (u,v) ] (3.6) 


Equations (3.4), (3.5) and (3.6) describe fully 


the transmission function of a Vander Lught filter. 


Before these equations are considered more closely 
it will be shown that this filter does indeed satisfy the 
requirements for a matched filter in that the impulse 
response of the filter contains the term s*(-x,-y) and this 
term can be separated from other terms in plane P3 cf 


Figure 2. 


Referring to Figure 2, when T(u,v) is introduced 
into plane P2, the inverse Fourier transform of T (u,v) 
appears in plane F3. Performing the inverse Fourier 
transformation on equaticn (3.4) (or equivalently Fourier 
transforming T(u,v) and then replacing x by -x and y by -y) 


shows that u(x,y), the light distribution in P3, is given 


by [17] 


u (x,y) =a26 (x,y) +S (X,Y) *S* (X,Y) 
tas! (~x,-y) *6 (x+w)f,y) tas (x,y) *6(x-w)f, y) (5.7) 


Convolution of a function with a delta function results in 
the same function with an origin shifted to the location of 
the delta function [18]. Equation (3.7) therefore contains 
s* (-x,-y) centered at x=-w)\f, y=0 and multiplied by the 


constant a. Furthermore, provided s(x,y) is limited in 
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Spatial extent this term can always, by an appropriate 
choice of w, be separated from the other terms of eguation 
(Sa/)i- The filter does therefore satisfy the requirements 


for a matched filter. 


From equations (3.4) to (3.7) some of the 


properties of the Vander Lught filter can be deduced: 


1) Comparing equations (3.4) and (3.7) it is seen 
that the first two terms in equation (3.4) do not ccntribute 
to the building up of s* (-x,-y), the important part of the 


impulse response. 


2) In order that T(u,v) remains non-negative a, 
the amplitude of the reference wave, must be chosen at least 


as high as the highest agplitude of S'‘' (u,v). 


3) The third term of equation (3.6) is both an 
amplitude and phase record of S(u,v) + S* (u,v). Here it is 
seen how the two complex waveforms are recorded together as 


a real record. 


4) If, during the recording of a filter, ¢s(x,y) is 
present in plane P1 of Figure 3 instead of s*(-x,-y) it is 


found that equation (3.7) changes to 


u (x, yY)=a26 (x,y) +S (x,y) *S* (X,Y) 
tas (x,y) *d(x+whf ,y) tas® (-x,-y) *6(x-wAf,y) (3.74a) 


The light distribution in plane P2 is therefore the same as 


before, except the locations of s(x,y) and S‘(-x,-y) have 
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been interchanged. This impulse response again satisfies 
the requirements for a matched filter and since s(x,y) is 
usually more readily available than s*(-x,-y), a Vander 
Lught filter is most often produced with s(x,y) in plane P1 


of Figure 3 instead of s'(-x,-y). 


3.2 The Lohman Filter 
A.W. Lohman and co-workers have published several 

papers about their method of using a computer and a 

computer-guided plotter to produce holograms which amongst 


other uses can serve as matched filters. 


The principle behind these holograms was first 
explained in terms of what the inventors have called the 
detour phase [8]. In a subsequent paper [9] a more 
detailed explanation of these filters was given and the 
function of the holograms examined in terms of diffraction 


theory. 


The description of the Lohman filter that follows 


is mainly based on the two above mentioned publications. 


The Lohman filter consists of a great number of 
rectangular apertures on an opaque background. The plane 
of the filter can be thought of as divided into a two- 
dimensional array of square cells with sides of length d. 
Each cell can be given two indices n and m to identify it. 
Let u and v be frequency coordinates in the plane of the 


filter and let the index n be increasing in the tu direction 
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and m be increasing in the #v direction. 


Referring to Figure 4 each cell contains one 
aperture whose width in the u direction is given by cd. 
Its height in the v direction is Wd. The aperture is 
centered in the cell in the v direction but displaced by 


Pima from the center of the cell in the u direction. 


Since the filter is being described in terms of 
frequency coordinates, all distances in the plane of the 


filter have a unit of 1/ (length). 


The variables (filter parameters) W,, and Pp are 
to be selected in such a manner that the mask has the 


properties required of a matched filter. 


The transmissicn function of such a mask can be 


written as [9] 


u- (Rte v-mnd 
T(u,v)= ) )Rect |——————— | Rect (3.8) 
n cd Wma 


) Y indicates summation cver all possible values of n and o 
nm 
and the Rect-function is defined by 


1 vibigal tho ag 69 Var) 
Rect (x) = (3.9) 
(th sie taper pips 
This filter is placed in plane P2 of the cptical 
system shown in Figure 5 and illuminated by the tilted plane 


wave derived from a point source located at X=X ) in plane 


Pl. 
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The problem at hand is to select the parameters 
Pom 2nd Wom if such a way that s*(-x,-y), the required 


impulse response, appears somewhere in plane P3. 


The system Shown in Figure 5 differs from that of 
Figure 2 only in that the point source in plane P1 is 
displaced off the optical axis in Figure 5 but is on the 
GxiSein rf .gutcr.. AS was pointed out in section 2.4 the 
displacement of the source off the axis results only in the 
complete light distribution in plane P3 being shifted. et 
therefore s' (-x,-y) appears in plane P3 of Figure 5 it does 
also appear in plane P3 cf Figure 2, only at a different 
location. Thus either system can be used to define a 


Matched filter. 


Assuming that the tilted plane wave has unit 
amplitude, the complex light amplitude of the wave emerging 


from P2 is egual to [9] 
B(u,v)=1 (u,v) Exp (j2™x,u) (3.10) 


The inverse Fourier transform of this amplitude appears in 
plane P3. This is given by 
u- (n+P 


DOR abot) Recta ee hect 
nm cd W 


Exp (j27x,u) 
(3.11) 


The above transformation is carried out in Appendix A and it 


is shown that b(x,y) is given by 
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b (x, y)=cd#Sinc[ (x+x,) cd] x 


L ) WymSine (Y Wynd) EXP {-J274[ (n+P,,) (x+x_) +my }} (3.12) 
where the Sinc function is defined by 


Sinc (x) =Sin (7x) 
TX (3.13) 
To establish that the Lohman mask is suitable for 
use as a matched filter, it must be shown that b(x,y) in 
equation (3.12), is at least approximately equal to 


s'(-x,-y) somewhere in plane P3 of Figure 5. 


Let it be assumed that s' (~x,-y) is limited in 
Spatial extent such that it is completely contained within a 
rectangular window w(x,y). The window w(x,y) has sides of 


length L, and L, in the x and y directions respectively and 


6 


is given by 

w(x, y)=Rect (x/L ,) Rect (y/L ,) (3.14) 
Then it is required to show that 

w(x,y) bd (x,y) =S' (~X-~Y) (3.15) 
The symbol = means “apprceximately equal to". 


In order to establish relationship (3.15) some 
approximations must be made to b(x,y) in equation (3.12). 


The approximations made by Lohman are the following: 
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1) Sinc { (x+x_) cd} =1 for Ix{SL,/2 
2) Sinc{i(yW,,d}21 for lytSLy/2, all hee (Seo) 


3) Exp (j2nxPd) =1 for [x{SL,/2, all eee 


Lohman has given both theoretical and experimental 
justifications for these approximations which are discussed 
in detail in reference [9]. Adopting these approximations, 
equation (3.12) becomes 


b(x,y)=cd2 } YW mEXE {-j2ma{ (nt+F) x tnx+my J} (SV7) 
n mM 


The above expression is seen to express b(x,y) in terms of a 
Fourier series, sc in order to facilitate comparison of 
b(x,y) and s"(-x,-y) the latter will first be expressed as a 


Fourier integral 
s' (-x,-y)=/{_S* (u,v) Exp {-j27(xutyv) } dudv (3.18) 


and finally s'(-x,-y) is approximated by a sampled version 


of the above infinite series 


s* (-x,-y) =d2 }) )S"(nd,md) Exp{-j2nd(xn+tym) } (3.19) 
nM 


The right hand sides of equations (3.17) and (3.19) can now 


be equated term by term which results in 
d2s* (nd,md)=cd@W Exp(-j2nd(nt+P yx] (3220) 


or if S* is written in amplitude and phase form 
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S*(nd,md)=A_ Exp {-j@_} (3.21) 


where Am and Dae are the amplitude and phase respectively 
of S*(u,v) at a point (nd,md) in. the u,v plane. Now 
equation (3.19) can be rewritten as 

etal MAS | (3.22) 


and 


Exp {-j@ J =Exp{-j27d (n*P_) x} (33723) 


This last equation can be simplified if XG is chcsen such 
that xd is an integer, Kk. Equation (3.22) can then be 
written as 


W A 


yen (3.24) 
Pom nm’ (27) 

These last two equations determine the values of 
oe te and asa such that the Lohman mask does act as a matched 
filter. It still remains to select values for the free 
constants K and c. The physical importance of these 
constants will be much mere easily understood when the 
relationship of Lohman and Vander Lught filters has been 
established. At this point it will only be noted that the 


values Lohman most frequently uses for these constants are 


K=1 and c=1/2. 


In summary, a Lohman filter is realized in the 
following way: 
1) The Fourier transform of the required impulse 


response is computed in a digital computer, usually 
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according to the fast Fourier transform algorithm due to 
Cooley and Tukey [19]. The value of the Fourier transform 
is computed for n x m points spaced on a regular grid. 

2) An array of n x m apertures is drawn on paper. 
The size and shapes of the apertures are described 
previously, in particular by Figure 4 and relationship 
(3.24). 


3) The plot is photoreduced to the required size. 


3.3 The Relationship he 


The Rela tween Van 


t der Lu 


er Lught and Lohman Filters 
It has been suggested [20] that the Lohman filter 
is the spatial equivalent of pulse modulation while the 
Vander Lught filter is based on amplitude and phase 
modulation of a Sinusoidal carrier. However, nowhere in 


the literature is it shown what the relationship between a 


Vander Lught and a Lohman filter is. 


In this section the transmission properties of the 
two filters will be analyzed and compared. The analysis 
will be conducted in the following way. The light 
amplitudes which expose a photographic plate during the 
recording of a Vander Lught filter are analyzed and S' (u,v), 
the part of the light amplitude which is derived from 
s'(-x,-y), is approximated by a piecewise continuous 
function (equation (3.28)). This allows the Vander Lught 
filter to be approximated by an array of cells, each having 
a Sinusoidally varying transmission function (equation 


(3.37)).- Finally these sinusoidal transmission functions 
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of the cells in the Vander Lught filter are approximated by 
binary transmission functions and they compared with the 
transmission functions of cells in a Lohman mask. The 
conclusion of the analysis is that the Lohman filter can be 
regarded aS a binary approximation to the Vander Lught 


LELter: 


In section 3.1 it was seen that the Vander Lught 
filter is the photographic recording of a complex amplitude 


given by 
C (u,v) =aExp (-j27wAfu) +S* (u,v) (3225) 


The variables are the same as defined in section 3.1 and 
S* (u,v) could alternatively be written as (see equation 


(325 ))) 
S* (u,v) =A (u,V) Exp {-j (u,v) } (3.26) 


Thrcughout the following analysis it must be born 
in mind that the spatial frequency coordinates (u,v) in the 
plane of the filter are related to the real coordinates 
(x,y) in that plane by u=x/Af and v=y/ 2. Thus the unit of 
the real coordinates (x,y) is that of length, whereas the 


unit of the frequency coordinates (u,v) is 1/ (length). 


Let it be assumed that s*' (-x,-y) is frequency- 
bandlimited so that S' (u,v) is zero outside some rectangle 
in the u,v plane and consequently the filter has finite 


dimensions. Thus it is assumed that the Vander Lught 
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filter is contained within a rectangular window E(u,v) in 
the u,v plane. This window is further divided intc nx an 
cells, each having sides of length d in the spatial 
frequency domain. Thus the window E(u,v) can be expressed 
as 


E(u,v)= . Ll perc Rect [yond] (3527) 


no d d 


Further, it is assumed that dis so small that S'(u,v) is 
almost constant within each cell, and the continuous 
function S*(u,v) is replaced in each cell by its value in 
the middle of the cell 


d 


S' (u,v) =) 2 S*(nd,ma) Rect fusna) Rect fea (3.28) 
nO 


With the above approximaticns the complex amplitude of the 
wave incident on the photographic plate in plane P2 of 
Figure 3 when the Vander Lught filter is being recorded, is 


written as 
C (u,v) =aExp {-j27wAf u} 


no d 


+) iS" (ndynd) Rect /ucad) Rect/ yond (3.29) 
d 


In equation (3.29) w is related to 6, the tilting 
angle of the reference wave indicated in Figure 3, by 
equation (3.2). In order that the reference wave has the 
same distribution across each cell in the Vander Lught 


filter, let @ be chosen in such a way that 
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wAfd=K (3.30) 


where K is an integer. 


By inspection of equation (3.29) it is readily 
seen that this choice of @ means that the periodic reference 
wave distribution has an integral number of periods inside 
each cell of width d. Furthermore, the absolute phase of 
the reference wave depends on the time-origin, which can be 
chosen at will. Let the time-origin be chosen in sucha 
way that the phase of the reference wave is ™ at the left 
hand side of any one cell. Since there is an integral 
number of periods inside each cell, this choice of time- 


origin ensures that the phase at the left hand side of each 


cell is si. Finally, a new coordinate x, is introduced for 
each cell. This new coordinate has an axis parallel to the 
u axis and x,=0 at the left hand side of each cell. Thus 


x is a coordinate commcn to all cells with the same walue 


of the index n but a different value for nm. 


In terms of the x, coordinates the distrikution of 
the reference wave can be written as 
aExp {-j2mwAfu}= )Rect fee sr act (3.31) 
n d 
Figure 6 shows the real fart of equation (3.31) plotted as a 
function of the coordinate x,. The plot extends over 
several cells in the u direction, and the value of K has 


been taken as 2. 
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Expressing the reference wave as in the right hand 
side of (3.31) allows it to be brought under the summation 


Sign in eguation (3.29) which results in 


C(u,v)=) )Rect fu-nd] Rect [v-md Gee (ay) (3.32) 
nM a a 
where 
Cum(Uev)=S* (nd,md) +aExp {-j27Kx /d-j 1} (3533) 


After development of the photographic plate the 


transmission function of the Vander Lught filter is given by 


T(u,v)= [C (u,v) [2 


) jRect ucnd) Rect [yoad] Can (Ue 2 (3.34) 
d a 


=_— 
—— 


n OM 


At any given point in the u,v plane, only one set of Rect- 
funtions has a non-zero value. Consequently, in squaring 
the sum of eguation (3.34) no cross terms appear. 
Furthermore, since a Rect-functicn only has the values zero 


and one, the square of a Rect~function is equal to the 


function itself. Thus expanding equation (3.34), T (u,v) 
becomes 
T(u,v)= ) Rect [urna] Rect /[¥-aa) [Cnn (OY) 2 (3.35) 
Dn O d 


Squaring the modulus of equation (3.33) yields 


Crm (Ue V) 2=a2taS* (nd,md) Exp{j2mKx,/dtjr} 


taS (nd,md) Exp {~j27Kx, /d- jr} (3.36) 
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If S*' is written in terms of amplitude and phase as in 
equation (3.21), T(u,v) can be written as 


d 


T(u,v)= ) )Rect a-n4) Rect was] 
D OD d 


UactAe #Zaa Cos {27 (Kx /dt1/2)-f 3] (Jaa) 


The transmissicn function expressed in eguation 
(3.37) is seen to vary in one direction only within each 
cell. This is the x, (cr u) direction. A cross-section 
of the transmission function indicated by equation (3.37) is 
plotted in Figure 7 as a function of a In ridure: 1.7K 


has been given the value 1 and ae is taken as 7/4. 


It was explained earlier that K determines the 
number of periods of the cosine-funtion within each cell and 
by inspection of (3.37) it is seen that 9g, determines the 
phase of the cosine-function. Referring to Figure 7, the 
distance Gumd by which the maximum of the cosine function 
deviates from the middle of the cell, and H,,,- the peak to 
peak amplitude of the sinusoidal transmission function, are 


given by 
G ane? nm/ (2 TK) and Hane 4aA nm (3.38) 


In secion 3.1 (c.f. equations (3.4) and (3.7)) it 


was seen that the constant term D given by 


D=a2+A2_, (S39) 
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did nct contribute to the part of the impulse response which 
was necessary for a matched fiiter. The importance of the 
constant D was to prevent the transmission function from 
becoming negative. If therefore the transmission function 
of each cell is reduced by a constant amount, such that it 
still remains positive or zero, the properties with respect 
to a matched filter have not changed. For this reason the 
modified transmissicn function shown in Figure 8 is equally 


valid as the one in Figure 7. 


Finally, in order to produce a binary filter it 
is required to approximate the continuous tone transmission 
function by a binary transmission function, i.e. an opaque 


mask with some transparent areas. 


There are many ways in which the sinusoidal 
transmission function of each cell can be approximated by 
one aperture in each cell. One of these is to use the 
arrangement previously described for a Lohman filter, and is 


drawn in Figure 4. 


If the arrangement of Figure 4 is used to 
approximate the transmission function of Figure 8, the 
following requirements should be met: 

1) The total amplitude of light coming from the 
aperture must be proportional to H_.- 
2) The aperture should be centered in the same 


location as the peak of the cosine wave, that is a distance 


G ard from the middle of the cell. 


add Yo soastzoqeh Sat. aee 
aotvonut acimaisaasit sit 970% 
st sede doe yfapoms taszenoo’ 6 ys 
toeqeo2 d3iw asitisqoie 3a7 aia i3 < 
od? pomde? allt 10% -bepasids rere prigrary: nee 
qifenpa ai 5 2aptt ot swore Ge 
¥ exept at diac as bitsy 


: teres iil é are ot sebi0 at (hban es 7 


sotoesnamgis eno! avowmitaod: eat ssanbiorin oi tela te 7 
Su pro ds .2.t «nolsoant eoteotaeabaz yout 5 yd wotgsan2 7 
chats snpasgewess somos Aviv tesa 7 


| : 
Istiowuate sds dokdw ai 2{ow yosm sms meee: {4 
ya ietaehscrags ad = liso é3se 36 noe sniet aozsslansoxs % 
eit ean 0% at sacdd Xo aa ».lies dons a& exusiays 0 a 7 

si bas \xes1/ apailod 6 10% Bsdizsesh yledoivany 2aseegneITe ji 


-» ommpet mk awezh 


oy boau 2t # sauplt Io soesepaiitey Sat 2r . 
aa? ,8 eames to notzoau7 ibtaabmeaszd adit _ 
item sd bivods stnensakupes W 
add nat painos sitptl ta ebextbir isjo2-20T a. 


esha od% nk bore235 5 antihadlentiadasie ss 
sonsialh # ek sod? .evew omteoo ait 30 Aseq oie e8 ai 


= 1 . oe 


a4 


Figure 8. 


TRANSMITTANCE 


4aA 
nm 


2aA 
nm 


Modified transmission function of a cell in the approximation to 


a Vander Lught? fitter 


45 


Since W,,, determines the total amplitude coming 
from the aperture in Figure 4 and Pim determines the 
displacement of the aperture from the middle of the cell in 


the u direction, these two requirements imply the following 


relationship between W,,, and Hy», and between P,,, and Gi 
Wom“Enm 
(3.40) 
Font Snm 


Since Hym%Anm 2nd Gym is equal to M,,/ (27K) it is seen that 
equation (3.40) expresses the same choice of filter 


parameters P and W,,, as Lohman arrived at by diffraction 


nm 


theory (equation (3.24)). 


In conclusion, the analysis of a Vander Lught 
filter has resulted in a kinary version of the filter which 
has the same filter parameters as Lohman derived by a 
completely different approach. Thus it has been 
established that the Lohman mask may be regarded as a binary 


approximation to a Vander Lught filter. 


3.4 Discussicn 

Now that the relationship between Lohman and 
Vander Lught filters has been established the free 
constants, K and c, of the Lohman filter can be examined in 


terms of a Vander Lught filter. 


As was indicated in section 3.3 the choice of K as 


an integer was to make sure that the reference wave has an 
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integral number of periods within each cell. This 
requirement leads to a simplification in the computation of 
the filter since the reference wave distribution is the same 
for each cell. K=1 implies that there is only one aperture 
in each cell which is obviously more convenient to plot than 
for example two apertures as in the case when K=2. 

However, as seen from equations (3.2) and (3.30) K is 
related to the tilting angle of the reference wave and thus 
also to the separation of s! (-x,-y), the required part of 
the impulse response, frcem the spurious parts (c.f. equation 
(3.7)). Thus K=1 is an appropriate choice only if the 


separation so obtained is sufficient. 


The choice of c=1/2 means that the cosine wave in 
Figure 8 is approximated by a square wave. This choice is 
intuitively easily understood. Other choices are, however, 


possible and have in some cases been made by Lohman [8]. 


A choice cf K=1 and c=1/2 leads to a problem which 
by Lohman is referred to as the overlap problen. To 
explain this, consider a cell with P,,=1/2. It is seen by 
considering Figure 4 that the aperture of this cell would 
extend into the adjacent cell by d/4. If now the adjacent 
cell happened to have P,,, , =-1/2, then the two apertures 
would completely overlap. This problem is easily 
understood in terms of the present analysis of the Vander 
Lught filter and can be corrected. Figure 9 shows the 


transmission function of a cell in a Vander Lught filter 
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which would lead to P,,=1/2, and it also shows how a Lohman 
filter approximates the transmission function. It is seen 
that the aperture extends into the adjacent cell. The 
problem of overlapping can be eliminated by choosing the 
alternative approximaticn to the Vander Lught transmission 
function drawn in Figure 9. Using this alternative 
approximation the aperture is divided into two parts, and is 


completely contained within one cell. 


It has now been established that the Lohman and 
Vander Lught filters are eguivalent, in that one can be 
regarded as a binary approximation of the other. The 
choice as to which filter is to be used in a matched 
filtering experiment does therefore depend primarily on the 
particular application and which of the two filters is more 


easily produced. 
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CHAPTER 4 


4.1 Introduction 

Despite the equivalence of a Lohman and a Vander 
Lught filter which has been explained in the previous 
chapter, there is an impcrtant difference in their 
production. While the Vander Lught filter is made using 
optical methcds only, the production of a Lohman filter 
requires that the signal to be matched, s(x,y), is sampled 
at regular intervals and the resulting array of numkters 


Fourier transformed using a digital computer. 


Prior to an algorithm due to Cooley and Tukey 
[19], known as the fast Fourier transform (FFT) algorithna, 
the Fourier transforming of large arrays of numbers was a 
very time consuming computation, so that the evaluation of 


Lohman masks would have teen highly unpractical. 


The FFT algorithm is a very efficient way of 
computing the Fourier transform of sampled data. A program 
written according to this algorithm requires as an input eT 
(n being an integer) equally spaced samples of a signal 
whose Fourier transform is required. The program then 
Fourier transformsthe array of samples, returning ay 
frequency components equally spaced in the frequency space. 
The highest (+ve and -ve) frequency obtained is 1/(2X) where 


X is the distance between samples, or in the case of time 
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series, X is the time interval between samples. 


Even when a fast Fourier transform algorithm is 
used, computing and plotting costs increase rapidly as the 
number of datapoints increases. The number of arithmetic 
operations required to ccmpute the Fourier transform using 
conventional techniques is proportional to N2 whereas using 
FFT they are proportional to 2Nlog,N [21]. N is the number 


of samples tc be transformed. 


Using the IBM 360/67 computer belonging to the 
Department of Computing Services, University of Alberta, it 
takes approximately 15 seconds (CPU-time) to compute the 
Fourier transform of a 128x128 array using the FFT 
algorithn. It would take approximately two hours to 
compute this transform using conventional techniques. For 
a 256x256 array the figures would be 1 minute and 36 hours 


respectively. 


Using the Calcomp 663 plotter belonging to 
Computing Services the plotting time for a 128x128 cell 
Lohman filter is approximately 1 hour, and for a filter of 


256x256 cells the plotting time would thus be about 4 hours. 


To give an indication of the cost of making a 
Lohman filter, the current rates for the computing 
facilities at this University are as follows: computing time 
(central processing unit) $600/hour, Calcomp plotter 


$16/hour. As a rough estimate, central processing unit 
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Storage costs during execution are similar to the cost of 


computing time used. 


The rates quoted above are the ones commercial 
non-university establishments are charged. At these rates 
the cost of computing and plotting a 128x128 cell Lchman 
filter can be estimated as $20 and the cost of a 256€x256 
cell filter would be about $85. The number of samples 
necessary to give a good representation of a signal is 
therefore an important factor in the production of a Lohman 


filter. 


In this chapter the sampling requirements for a 
Lohman filter will be examined, using fingerprints as test 


patterns. 


The reason for using fingerprints as test patterns 
is, as menticned in Chapter 1, that the work done on 
fingerprint identificaticn using Vander Lught filters [5,6] 
indicates that fingerprints are well suited for detection by 
matched filters and unlike manmade objects, such as 
alphanumeric characters, fingerprints seem to possess the 


arbitrariness ccmmcn to many naturally occurring objects. 


4.2 Sampling Theorem = Aliasing 

It is required to use the fast Fourier transform 
algorithm to compute the Fourier transform of S(x,y). To 
this end s(x,y) is sampled over a regular grid of sampling 


points, spaced a distance X apart in the x direction and a 
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distance Y apart in the y direction. The array of samples 
a(x,y) so obtained can be written as 

a(xey)= ) )s (x,y) 6(x-nX,y-mY) (4.1) 

n m 
The FFT algorithm can be used to compute the 

Fourier transform of a(x,y) which is also known as the 
discrete (or finite) Fourier transform of s(x,y). The 
question arises, how accurately the discrete Fourier 
transform approximates the continuous Fourier transform of 
S(X,Y)- The answer to this question is given by sampling 


theory: 


If s(x,y) is a continuous function and a(x,y) is 
related to s(x,y) by equation (4.1), then the frequency 
spectra A(u,v) and S(u,v) of the two functions are related 
by [18] 


A(u,v)=[1/(XY¥) ] } )}S(u-n/X,v-m/Y) (4.2) 
n MQ 


The right hand side of equation (4.2) is propcrtional to an 
array of functions equal to S(u,v) centered at (n/X, m/Y) 
where n and @ can have any values between -~ and © Thus 
the nearest S(u,v) functions are separated by 1/X and 1/Y in 


the u and v directions respectively. 


If s(x,y) is a bandlimited function then S(u,v) is 
non-zero only in a limited region of the frequency plane. 
Let it be assumed that s(u,v) is non-zero only for 


frequencies j{uj<W and {v{<W. Under these conditions the 
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various S(u,v) functions in equation (4.2) are separated 


provided X and Y are chosen such that 
X<2W and Y<2wW (4.3) 


The situation described above is drawn in 
Figure 10a for the part cf the spectrum along the line v=0. 
Only the central part of the infinite array of S (u,v) 


functions is shcwn. 


Suppose now that X is increased such that the 
relationship (4.3) is no longer true. Then as is shown in 
Figure 10b the various S(u,v) terms are no longer separated 
but overlap. Thus it is no longer possible tc recover 
S(u,v) exactly from A(u,v). The damage to S(u,v) is 
referred to as aliasing and it is said that high frequency 
components from an adjacent S(u,v) function alias high 


frequency components of the central S(u,v) function. 


Pricr to the experiments with fingerprint 
identification described here, some experiments were done 
with Lohman filters in character recognition. In these 


experiments the problem of aliasing was quite apparent. 


Figure 11 shows the amplitude of the Fourier 
transform of the letter A, computed using the FFT algorithm. 
The computed Fourier transform is composed of 128x128 
datapoints. Each of the datapoints is represented in 
Figure 11 by a black square, the area of the square being 


proportional to the amplitude of the Fourier transform at 
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that point. 


When the computed Fourier transform of Figure 11 
was compared with the optically obtained Fourier transform 
of the same letter A, it was found that the line of spectral 
components indicated by (1) in Figure 11 was in fact the 
continuation of the line (2) from an adjacent spectrum, 


appearing in the central spectrum because of aliasing. 


Intuitively one might suspect that fingerprints, 
being complicated patterns, possessed a wide range cf 
Spatial frequencies. If this was the case, it would be 
necessary to employ a great number of samples to adequately 
discribe a given area cf a fingerprint. This would lead to 
high computing costs and thus render the Lohman filter very 
uneconomical for fingerprint identification. For this 
reason a survey was undertaken to investigate if 
fingerprints are indeed bandlimited, and if so, to establish 


an appropriate sampling rate. 


In order to investigate the frequency content of 
fingerprints, ten fingerprints were selected at random from 
a file of prints, originally obtained from the R.C.M.P. 
Each of these prints was magnified approximately 10 times 
and the resulting photograph sampled at a high rate to 
minimize aliasing. The resulting array of samples was 
Fourier transformed by the FFT algorithm and the resulting 


frequency spectrum analyzed. 
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The sampling of the enlarged print was performed 
using an image conversion system belonging to the department 
of Computing Science. This system consists of a Sony DXC- 
2000A TV camera which scans the print, and a quantizer 
designed and built by Technical Services, University of 
Alberta. The quantizer samples the video output of the TV 
camera at instants that correspond to a regular grid of 
points on the fingerprints. The quantizer is connected to 
a PDP-9 computer which stores the samples from the quantizer 
on a magnetic disc. Also connected to the PDP-9 is a 
cathode ray tube display unit which allows a picture, stored 
in a sampled form on the disc, to be redisplayed and 


inspected. 


The guantizer is capable of distinguishing between 
64 levels of gray and the totai field of view of the TV 
camera corresponds to 44€x448 samling points. Smaller 
areas can be selected with a corresponding reduction in the 


number of sampling points. 


The time taken by the system to sample and store 
one picture is about 40 seconds. Thus the system is well 


suited for the processing of still pictures. 


The samples of the prints were very closely spaced 
or over 20 samples per linear millimeter of the actual 
fingerprint (over 400 samples/mmé?). This very high rate of 


sampling was chosen in order to minimize aliasing so that it 


to _tiewaviaD ,aeolvzst ieoinioet ss siahs elaine 
Vv? ogf Yo tugtle bebiv sf? colguae Tesetnsup of ad aLN 

to bite asivesa 2 o4 baneeetson 08? atopseet*se'ieebeo. - 
oF batoynnod at absitnivy sil . -aembagzepag? oes 0 aah. | 7 
19nL*daWp sd¥ api aaiquéa edd aoante dotde wasugnon 0g a : - 


2909 sit of be taenoY oala paib eee . 


betade .oaéFois « avaltes dabée fiav yelqeth edus yes sbostea — 
yn hevelgeiber ad py ,9ath eA? a0 ated be Iqase sab 


. ‘ 
: >” 


-_ 


>a ese parhuecvpatseth fo sidajay 25 tezivnsep s4T bs Saas : 

V? af? Jo weilv to bier! intod Bdd bas ye7p To elsvel ¥o 

teiieuz  .atiboq pwiiase 8ode9 80 Of ahdoqgesatod ‘ezedso. : 

ag3 Wi aokiouber onibuege orien i7 iv petoeles ed aso Seats _ 
+FIgiog en tiques 4 re 

e10t2 fds alygwse oF aetaye ode Yd aatet saso GAT 7 = 


[lov ct abwaye okt avd? .,abwouwe 08 goods al esessigueto. — | 
_ desu ty Live to ¢uleseoony o¢2 402 Beekee | 

beoeqe piezolo yasv s14¥ etoksg sds Io aaiquee aT we. t 7 
SAUT26 REE Ae Teh, SRAM, AES See RE | a ee - 
aerorien eat ae! see paatda es (000, seep pe 


11 dete om cotaalts sxiutacs 9s seb20 ab aomMe au 


58 


could be neglected. 


The array of samples obtained in this manner was 
then Fourier transformed using the FFT algorithm resulting 


in a frequency spectrum cf each of the ten prints. 


AS was just mentioned, it was assumed that 
aliasing was negligible. This assumption was justified by 
the fact that the highest spectral components obtained were 
very small compared to the lower ones. The computed 
amplitude of these high ccmponents is the sum of the 
component in the fingerprint at that frequency and the 
amplitude which appears at that frequency because of 
aliasing. Since the total amplitude computed for high 
frequencies was low, certainly the amplitude appearing at 


these frequencies because of aliasing must have been low. 


In order to give some representation of the 
bandlimitaticn of the prints, the total power contained 
within the frequency band f{u|<W and (v{<W was calculated for 
several values of W. All 10 prints showed similar results, 


so only the two extremes are presented here. 


Figures 12 and 13 show the two extreme cases. 
The print represented in Figure 12 is the print with the 
greatest concentration of power at low frequencies (i.e. the 
line plotted in Figures 12 and 13 rose fastest for this 
print) whereas the print represented in Figure 13 is the 


print out of the ten with the least concentration of power 
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Percentage of total power as a MincwLont of 
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Low-frequency print. 
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at low frequencies (the line in Figures 12 and 13 rese 
Slowest for this print). From now on the former will be 
referred to as "the low-frequency print" and the latter as 
“the high-frequency print". From Figures 12 and 13 it is 
seen that the prints are indeed fairly well bandlimited. 
Even in the case of the high-frequency print, 90% of the 
total power of the print is contained at spatial frequencies 


lower than 4.2 lines/mn. 


In selecting an appropriate sampling rate it 
should be kerft in mind that the power at frequencies outside 
the band-limit W dictated by the sampling rate 
(c.f. equation (4,3)) will net only be neglected, but will 
in fact, through aliasing corrupt the frequency components 
within W. On the other hand it is desirable to keep the 
Sampling rate down so that a large area of a fingerprint can 


be represented by a small number of samples. 


With these criteria in mind it was decided that a 
sampling rate of 11 linesy/mm (yielding a frequency hband- 
limit of 5.5 lines/mm) was an appropriate compromise. 

Using this sampling rate, in the worst case (see Figure 13) 
only about 5% of the total power in a print is affected by 
aliasing and in the best case (see Figure 12) about 2.5% of 
the total power is affected. At this sampling rate 
1.2x1.2 cm2 of a fingerprint can be represented by 128x128 
samples, and 256x256 samples would be sufficient to 


represent a complete fingerprint. 
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PRODUCTION OF A LOHM 
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5-1 Data Collection 

in Chapter 3 the theoretical basis of the Lohman 
filter was explained and in the last chapter it was 
described how a pattern can be sampled in order to obtain 
the data required to construct the filter, In this chapter 


the complete process of making a Lohman filter will be 


discussed. 


Figure 14 shows a flow diagram of the operations 
required to froduce a Lohman filter, starting with a 
photographic record of a fingerprint and finishing with a 


Lohman filter matched to that print. 


The fingerprints used in the experiments were 
obtained as records on a microfilm. Thus the first step in 
the producticn cf a Lohman filter is to enlarge the print 
whose filter is required onto photographic paper. The 
scale of the repreduction is determined by the image 
conversion system used tc sample the prints and was in the 
present work approximately 10 times the natural size. This 
enlargement was performed by Photographic Services, 


University of Alberta, and was done by a two-step process. 


The image conversion system used was described in 


the previous chargter and is drawn in a schematic form in 
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Figure 14, The image ccnversion system outputs an array 
humbers describing the darkness of the print at points on 
regular grid. In the experiments 128x128 samples were us 
and a sampling density of 5.5 samples/mm in the two 


perpendicular directions in the print plane was used. 


The output used from the PDP-9 computer was in t 
form of punched paper tape and contained the data asa 
128x128 (=16384) element vector of integers 0 to 63, 
corresponding to the 64 levels of gray distinguishakle by 


the quantizer. 


Unfortunately the computing system at the 
Department of Computing Services did not include a high 
speed paper tape reader at the time of the experiments. 
order to input the data to the IBM 360/67 computer the pap 
tapes first had to be translated onto magnetic tape which 
could be read by the IBM 360/67. This translaticn was do 
in the Divisicn of Biomedical Engineering, Clinical Scienc 
Building, University of Alberta, using their ccmfuting 
facilities. The translation is not indicated in the flow 


diagram in Figure 14. 


5.2 Computing and Plotting 
Once the array of samples has been stored in the 
IBM system it is required to evaluate the parameters 


describing a Lohman filter and then use the Calcomp 663 


plotter to plot an enlarged version of the filter. 
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Twc computer programs were written for this 
purpose. The first program uses the Fourier transforming 
Subroutine HARM which is available in the IBM Scientific 
Subroutine Package to Fourier transform the array of samples 
obtained from the image conversion systen. The elements of 
the Fcurier transform are then converted to an amplitude and 
phase notation instead of the real and imaginary notation 
outputted from HARM. Finally, the amplitudes and phases 
are used to calculate the parameters of the filter according 
to eguations (3.23) and (3.24). Both this program and the 
program used to create the plots are listed and explained in 


Appendix B. 


As pointed out in Chapter 3 there are many Shapes 
of apertures that fulfill the requirements placed ona 
Lohman filter. Referring to the notation of Figure 4 the 
shape of apertures used in the present filters was the same 
as drawn in Figure 4 except that W,, was not varied but was 
taken as a constant equal to 1. Instead c was allcwed to 
vary ketween 0 and 1/2 tc determine the amplitude of light 
emerging from each cell. Thus each aperture extends over 
the whole ceil in the v direction and the width in the u 
direction was varied from 0 to 1/2 the cell width. If an 
amplitude calied for an aperture width of less than one 
stroke of the plotter pen, a single line was drawn and the 
length reduced appropriately. If the area of the required 


aperture was less than the pen tip no aperture was drawn. 


iofunes ts glee aaenpontaiene ae 
to omessis adv 6. eeseye noisrovass syake-sar aoxt benkaddo. . 
hie shyrelqas as oF hatteraop asd? ots szoteaes? 13104 ods ‘2. 
netteton Yindipset bas lesa edt te Bemnins oat | tt . 
seckdg Bos 2atwttlqms edt ,ylisors NEAR O22 sessagsvo 
entbrouss 1ssfbi, edd Yo 2usteubisg Sat sinteein ot bene oe 


ot? bas men bdr tog = (PSE) Bas (6S.€) anoeere iy : 


ag vane ‘pas betel oma ajoly sit atse7> ot beas sstpomq 
/ | +«@ xibaegga 


4 
eeyjkhde ¥itee Ste stsa? E Ved ed9 al IVG Bstatog aa 7 


“wit 
s no Beostly Stpemethupes. sat [let ivr sede eniss30qn 26 
si) § azupi) Jo soltston sd3 of Pa lrIeien ~nerizt asadot 
onec off Hey aintili tasecrg sds gi hoso Some seqs to ogece 
eee S04 hodtiny zon an, 7802 ?4e0ts 6 saepla “ni awezb ae 
oS Povolls -ow o bastent  .f of Lavpe twagtaos 4 as nodes wed 
tépel to shutilyes sd? Sathtetsb 29 SNF one 0 neowsed en 
aavo gaerxe siudtsqgs dnas eud? Jiisp doze soz? pelesane 
uv edt al dtbiv odt ban doLd9szi6 v odd ai IDS0 elodw ods 
as it .itbiw Liem eu? S\f of @ mozt-befesy eaw soksoeaeh 
$9 gett Zeal lo Withtv saudidgs oe wOX boften hot kigas 
odd bas aueah cav oAkl sibake © .09q aeftoly bat 30 stows . 
hoyiupes ett lo Kerr alt Mt | .ylereksqgongge Bevetsd Aspaed 
ANSTO GS¥ SIUsIsds On YLT Mey 212 nade Geel gow eapPreqe | ; 


idl 


66 


AS an indicaticn that neglecting these small 
apertures was reasonable, it should be pointed out that the 
largest aperture neglected is only about 1/500th of the area 


of the largest aperture. 


The plotting program listed and explained in 
Appendix B picts the filter column by column, decides for 
each cell which type of aperture is required and draws it 
according to the procedure described above. This 
plotprogram makes use of subroutine PLOT which is available 


in the Calcomp Subroutine Package. 


The final step in the production of the filter is 
to photographically reduce the plot onto photographic film 
such that apertures drawn in black on the plot become 


transparent and the white background becomes opaque. 


The scale used in the plotting was .22 inches/cell 
(.56 cm/cell) and this was reduced to a filter of size 100 
microns/cell or a 55 times reduction. This reduction was 
again carried out by Photographic Services and was done in 


two steps. 


5.3 Cost of a Filter 

In the last chapter the cost of computing and 
plotting was found to be about $20 for a 128x128 cell filter 
and about $85 for a 256x256 cell filter. Costs other than 


computing and plotting costs are the cost of sampling a 


picture, and that of the photographic processing. Tteis 
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difficult to assess the cost of sampling since in the 
present work the image ccnversion system was lent free of 


Charge by the Department cf Computing Science. 


The ccst of the photographic work involved, both 
the enlargement of the print and photoreduction of the 


plotted filter, was about $5-10 for each filter. 


Thus the cost of a complete filter can be 
estimated as $25-30 for a 128x128 cell filter and about $95- 


100 for a 25€x256 cell filter. 
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CHAPTER 6 
EXPERIMENTS 


6.1 Experimental Arrangement 

In Chapter 2 spatial filtering is explained in 
terms of the optical system shown in Figure 1. This systen 
can be used for any spatial filtering experiment and has the 
advantage that Fourier transforming properties of lenses can 
easily be explained in terms of this systen. The system 
sketched in Figure 15 is, however, better suited for 
practical experimentation [23] as will be seen later and 
will be used in the experiments described in this chapter. 
The properties of this latter system are derived in detail 


in reference [23] and only the most important results will 


be pmentioned here. 


In order to obtain an exact Fourier transform of 
any object it must be placed one focal length in frent of 
the Fourier transforming lens and then the exact Fourier 
transform appears in the back focal plane of the lens [14]. 
This is the arrangement used in the system shown in 


Figure 1. 


If on the other hand, the object is placed behind 
the lens less than one fccal length away, the Fourier 
transform of the object still appears in the back fccal 
plane of the lens but this time it is multiplied by a 


spherical phase factor [ 14]. The filtering system drawn in 
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Figure 15 makes use of this arrangement of lens and object 


to perform the Fourier transformation. 


Referring to Figure 15, L1 is the lens which 
performs the Fourier transformation and is illuminated by a 
plane wave. The object to be filtered is placed in plane 
P1 which may be moved alcng the optical axis of the system. 
AS previously mentioned in this case the Fourier transform 
of the object which appears in plane P2, the back focal 


plane of lens L1, is multiplied by a spherical phase factor. 


Using the notation of Chapter 2 and with reference 
to Figure 15, the spatial frequency coordinates, u and v in 


plane P2, are related to the real coordinates by 


u=x/ (AD) and v=y/ (AD) (6.1) 


Comparing the above equation to equation (2.17) it 
is seen that they are the same except the focal length f in 
equation (2.17) is replaced by D, the distance between P1 
and P2 in Figure 15. Since plane P1 can be moved along the 
optical axis, D is a variable and thus the scale of the 
transform in plane P2 is variable. This feature is one of 
the main advantages of the system shown in Figure 15 over 
the one in Figure 1 since it allows errors in scale of 
either the object or the filter to be compensated fcr by an 


appropriate choice cf D. 


Ancther advantage of the system in Figure 15 is 


that it is more compact since the Fourier transformation 
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PLANE 
WAVE Le 2 y 


Figure 15. Optical filtering system used for experimentation. 
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only requires cne focal length compared to two focal lengths 


required by the system of Figure 1. 


Referring again to Figure 15, lens L2 images the 
filtered version of the pattern in plane P1 onto plane P3. 
In order that a focused image appears in plane P3 the 
relative positicns of P1, L2 and P3 have to satisfy the well 
known imaging condition [14] which in terms of the 


parameters in Figure 15 can be written as 
fa + 1sb=1/f5 (6-2) 
where f, is the focal length of lens L2. 


6.2 


la) 


ositioning Accuracy Required 


It can be seen by inspection of equation (6.1) 
that the origin of the spatial frequency coordinates in 
plane P2 is defined by the optical axis and does not depend 
on the exact position of the object perpendicular tc the 
optical axis in plane P1. The positioning of the object in 
this direction is therefore not very critical. It is only 
necessary that the object remains fully illuminated during 
the filtering experiment. However, as the object is moved 
in plane P1 (perpendicular to the optical axis), the 
filtered image also moves in plane P3. Thus in order to 
decide whether recogniticn has taken place or not the point 
of light indicating recognition must be searched for in 


plane P3. 


It is, however, important that the filter is 


a4t expewi SJ anel et 
6% suatg dido fa weekg. Cf 
sd+ £9 sesfq at Pee F Beas 
Iiew ed? yRedsee 09 avad £4 Bre, Sa pipf avg ee 


odd Yo awsse nf dod (OC), mobtbinos eubeont 
as asttiay sd oko éf oamplt of 


($8) Zo | 2\Psaiyt + Ne 


ve Si anel 20 ddemet Ieont sd? at a ose 


‘bealuees seemed paganbttand 6. 


(f .4) nosatna to nogtosgent yd asec sd a6D 22 | 

ak 233snlbitoos yonsupes7t Lesseqa edt 30 aigizo odd suas 
ineqeh ton asch Sun Sixes levidgo sda yd /beltieh eS 24 

att oF 1sivstbhasqaay toutdo ad? 2o aotsleoq toste ao 
nz Joetdo sd¢ to BeLadisteog eat .f9 saslg ui alee isoidqo | 


ae a 
a a sk lr 


ee 


yino 2% #3. sievursty> yiev Feit eagdensils at solsoenth - 
gaiavh boteniwnlit yilvui seiouer tostdo Sdo cede ; 
hovon ef tyatdo edd as ,seveWoR .tnewi19qes owe tes£b2 edt : 

ait \(aixn fsottqo eds o2 - s9.1psbneusegh ra onekg ah *) 


: 
0? ptt gi Bud? -E4 oaelg ak 2¢vou oats aggnd Beuediit ia 


tetog edo Jom 20 aonly aad59 Sei gotTiNpoos. aedssde sbioeb . 
fs sot tofassenis Selene se RO ane 
Mains 


‘ 


et aea{th odd 267 tneraoges .aevewOd yeh aT 


W2 


accurately placed in plane P2. As previously mentioned the 
origin of the spatial frequency coordinates is the optical 
axis. In order that the filter matches the Fourier 
transform in plane P2, the point on the filter which 
corresponds to zero frequency on the Fourier transferm of 
the object must be placed on the optical axis with 


sufficient accuracy. 


In order to obtain some estimate of the accuracy 
required in the positioning of the filter perpendicular to 


the optical axis the following argument is adopted: 


In the derivation of a Lohman filter it was 
assumed that the complex transmittance varies slowly within 
each cell. (This was the assumption which allowed the 
transmittance of each cell te be approximated by one 
aperture only.) if therefore the error in placement of the 
filter in plane P2 is small compared to the dimension of a 
cell, the resulting error in light transmitted through the 


filter can be expected tc be small. 


In the experiments performed all the filters 
consisted of 128x128 cells and had an overall size of 
13x13 mm2. Thus the side of each cell is approximately 


0.1 mm or 100 microns. 


Let the phrase "small compared to the dimension of 
a cell" mean “less than 1/10 the cell side". Then it can 


be estimated that it is necessary for recognition that the 
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filter deviates no more than 10 microns from its correct 


position. 


Of course it is possible that the complex 
transmittance does not only vary slowly within each cell, it 
may vary Slowly within several cells. If this was the case 
the above estimate of the accuracy to be achieved is a worst 


case estimate. 


Vander Lught has published some theoretical 
studies about the deterioration in signal to noise ratio as 
a function of the error in the placement of his filter [24]. 
These studies show that displacement of the filter in the 
axial direction is relatively unimportant so that the 
accuracy required can easily be achieved by manually sliding 


the filter hclder along an optical bench. 


In designing a holder for the filter the main 
criterion to be met was therefore that the filter could be 
positioned with an accuracy of at least 10 microns in the 


directions perpendicular to the optical axis. 


Referring once more to Figure 15, the plane wave 
used as illumination was derived from a Spectra Physics 
Stabilite Helium-Necn laser of wavelength .6328 micron and a 
Spectral Physics model 311 collimating lens assembly. The 
Fourier transforming lens L1 had a focal length of 152 cm 


and a diameter of 10.4 cm. The imaging lense L2 hada 
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focal length of 54 cm and a diameter 5.2 cm. The two 
lenses L1 and L2 were mounted on a heavy cast-iron optical 
bench together with a holder in plane P1 containing the 
object (fingerprint) and a holder in plane P2 for the 
filter. Ail the components on the optical bench could be 
slid along the bench and by means of a vernier positioned 


with a 0.2 mm accuracy. 


Plane P3 contained a simple paper sheet to enable 
the inspection of the light distribution in the plane. All 
power measurements in plane P3 were performed by the use of 
a Spectra Physics model 401B power meter. In all cases the 
power was measured through a circular aperture of diameter 


approximately 1 mn. 


The filter holder was made at the Machine Shop, 
Department of Electrical Engineering. Figure 16 
illustrates the mechanical design used to achieve the 
required accuracy in positioning the filter. The outside 
frame is composed of two frames . Each frame is cut out of 
3/4 inch sheet aluminium and groves are cut into the frames 
to accommodate the hardened ground steel rods shown in 
Figure 16. The sledge shown slides along the steel rods on 
precision low friction ball bushings arranged in the fashion 
indicated in Figure 16. In order to control the position 
of the sledge accurately a micrometer is mounted through the 
outside frame tc push the sledge along. A steel spring 


holds the sledge against the micrometer. 


ees oft ven S62 
Lsstsyo rorinaaea yond 6 a0 besauga $3 


eit entutssaco 1 smald ah 2s peg «al aaas | 
tiaras 


eit sot £4 suaig ab 2 
Lgand @ i 


cv ie 
a 


ed bives domed Isottqo oft 20. 
ome han ae 7 


sideas oF teeda iqnq shyaie » betkstaop £9 saat 

fia ,eeglg sdi/nt apitodiatetb Idpkl sdt Yo aoisoeqeat odd 
to sau 949  begrottsg ote €9 sueig at eteae wasee 2ovog 
Sit 9B 82 the at ,Isfou voy 2900 Ishoe satay dd sas0qe 5 
panei en sinttteqe isfuotto s dpaozd? Satpaces asw aun 
«om | ¢fossaixozqgs 

,qod2 snbfoen eae te shat esw aehlod asthit sdk” 
af saul .patieeaivak Lsoitsoalad Yo teens z6qGed 
sit susidoe oF beau. apezeb isvinsdzea of? ee2zsitse0Lit 
shietiio eit .retiti sd? ouinoitteoy ni YosI99os poxtupea 
to tvo go> @t exert Wont 8 « semest on? Zo Beeogaos be oasz3 —— 
souett odd ofat fun 925 asvozt bas solaiwils Joode ‘teak axe 
nt awode abot Iseta bavosp Beaebrsd od? ot sicmee226 of 
ao 2bor Ieete eit pools sabile nvode epbele sit -8F sipplt 
noldck? out ak bOpasT2S epaidoud Ifed noises? wol nobetoesq 
ooktiaog sat Lograop of ashzo et «dF etmplt at betsotbat 
oft dpvotds baynuem ef 2otesor0Le 6 yless2n908 spheke 94% to 
vaizqa Ieote A -phols spbels od? dang of amez2 ebletoo 
jzerenorsia of3 Santspe epbese edt mblod 


ie 


*uOTRJOSATpP N UT ASAITTF uot#tsod 03 ubtsep TeoTueYyooW *9T oanbhty 


dod 


ONTHSNd TIVd 


ONTads 
aOdgad1S 


-nofigoe tid vw nt yexliG noktiecq of 


npiesh levinedost 


76 


The micrometer is a 40 turns/inch micrometer which 
indicates that each turn cf the micrometer shaft corresponds 
to a motion of the sledge of approximately 625 micrens and 
ten microns ccrrespond tc 6 degrees of turning which is 


easily controlled. 


Finally, in order to achieve the same accuracy of 
positioning in the cther direction a similar device is 
mounted on the sledge, the orientation of the seccnd device 
being perpendicular tc the first one. The apertures on the 
two devices overlap to eS ieweniahe to pass through. Now 
the filter can be mounted on the sledge of the second 
device, covering the circular openings of both sledges. 

The filter can be mounted with or without a liquid gate 


Which will be described shortly. 


When a Vander a Lught or a Lohman filter is said 
to affect only the amplitude of a light wave it is assumed 
that the photcgraphic film used to record the filter is of 
uniform thickness. With commercially available films this 
is usually nct the case and the variation in thickness of 
the film can give rise tc errors in the phase distribution 
of a light wave emerging from the filter and reduce its 
performance [17]. In order to reduce the effect of 
thickness variaticn of the film the previously mentioned 


liquid gates are used. 


A liquid gate consists of two highly polished 
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glass plates between which the film is inserted together 
with a small amount of oil with index of refraction as close 
as possible to that of the film. In this arrangement the 
oil fills the depressions in the film and the spatial 


variation in optical path length through the gate is small. 


The glass plates used for the liguid gates were 
polished in the optics shop of this Department and dibutyl 


phtalate was used as index matching oil 


AS previously stated, liquid gates could he 
mounted on the filter holder and also on the holder which 
contained the object (fingerprint). The latter holder 
could also be tilted about the optical axis and the angular 
position set accurate to +10 minutes. This enables the 
object to be adjusted in angular position to match the 


Pricer. 


6.4 Experimental Results 
The first experiment was designed to confirm and 
demonstrate the ability cf the Lohman filter to recognize 


fingerprints. 


Eight prints were selected from the file cf prints 
acquired from the R.C.M.F. of these eight prints three 


were used as originals tc produce Lohman filters. 


Each of these 2 filters was then used to detect 


the appropriate print frem the selection of eight. 
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The prints were prepared in a form suitable for 
the experimental system described in the last section. A 
1.2x1.2 cm2 section of the print was made into a negative 
transparency of 0.8x0.8 cmé2. This particular reduction 
(33%) was chosen because it was convenient to work with in 
the system used. These transparencies were used as objects 


in the experiments and mcunted in plane P1. 


The filters were made by magnifying the 
appropriate object about 15 times (10 times the real size of 
a fingerprint) cr to a size suitable for sampling by the 
image conversion system described in section 4.3. The 
128x128 array of samples was Fourier transformed anda 
Lohman filter made using the programs and computing 


facilities described in Chapter 5. 


Each of the 3 filters was then tested by the 


following method: 


Both the filter and the object which it was known 
to match were mounted in liquid gates and inserted into 
their holders. The light intensity distribution in plane 
P3 was observed and the position of the filter adjusted by 
use of the micrometers on the filter holder, and by sliding 
the holder along the optical bench, until a recognition spot 
was found. The system was then further adjusted so as to 
maximize the power in the recognition spot. After these 
adjustments the system is kncwn to be correctly set for 


identification of prints prepared in the fashion described 
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earlier. (The system is set for the correct scale of 
objects and the zero frequency in the filter is at the 


optical axis.) 


Now the print in plane P1 was replaced by cne of 
the other eight and a recognition spot visually searched for 
in plane P3. If a recognition spot was found its fower was 
measured but ctherwise the power in plane P3 where the 
recognition spot should have appeared was measured. This 
was repeated for all the eight prints, ending with the one 


the filter was known to match. 


The results of these experiments are shown in 
tables 1,2 and 3. The tables indicate for each filter and 
each object whether a recognition spot was observed in plane 
P3 and the power in plane P3 at the location where the 
recognition spot should have appeared. In each table the 
power has been normalized to the power of the recognition 


spot. 


From the tables it is seen that a recognition spot 
is visually cbserved only when the object is the print to 
which the fiiter is matched. The maximum power in the 
location of the recognition spot when the filter does not 
match is only 14% of the power measured when matching to the 


filter occurs. 


Next the effect cf a small displacement of the 


filter in a direction normal to the optical axis was 
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RECOGNITION POWER 


OBJECT SPOOL IN 
NO. OBSERVED PLANE P3 
1 YES Poo 
2 NO Shad Oe 
3 NO 0.04 
4 NO OU 
5 NO Os O es: 
6 NO 0.04 
7 NO pais @) 
8 NO Ou 
TABLE 2. RESULTS ‘OF (FINGERPRINT ‘RECOGNITION 
EXPERIMENT. FILTER MATCHED TO 
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TABLE 2. 


RECOGNITION POWER 
OBJECT SPOT IN 
NO. OBSERVED PLANE 

a} NO 0.14 

2 YES 1.00 

3 NO 0.04 

4 NO 0.06 

5 NO 0.06 

6 NO 0.02 

7 NO 0.06 

8 NO 0.04 


RESULTS OF FINGERPRINT RECOGNITION 
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TABI 3. 


RECOGNITION 
OBJECT pion 
NO. OBSERVED 

£ NO 

2 NO 

3 LES 

4 NO 

5 NO 

6 NO 

of NO 

8 NO 
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OBJECT NO. 


POWER 
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PLANE P3 
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studied. The filters were first adjusted to maximize the 
power in the recognition spot and then deliberately moved 
from their correct position in both the positive and 
negative x and y directicns. The reduction in power in the 
recognition spot was reccrded as a function of the 


displacement. Figure 17 shows a typical result. 


From Figure 17 it is seen that the power in the 
recognition spot drops by 10% if the position of the filter 
is incorrect by approximately +75 microns. Thus the 
estimates of the accuracy required in positioning presented 
in section 6.2 are seen to be too stringent. This could 
have been expected since the estimate applied to the worst 


case. 


The effect of angular displacement of the filter 
from its ccrrect position was also studied. In these 
experiments the object holder was used to initially adjust 
the relative angular position of print and filter to 
maximize the power in the recognition spot. Then the 
angular position of the frint was changed in both the 
clockwise and counterclockwise directions and the power in 
the recognition spot recorded as a function of the angular 


position. 


Figure 18 shows a typical set of these 
measurements. If a power loss of 10% is once more taken as 
a criterion, it is seen from Figure 18 that the accuracy of 


positioning in the angular direction has to be better than 
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RELATIVE 
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Figure 17. Power in the recognition spot as a 
function of placement error perpendicular to 


optical axis. 
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approximately +1.5 degrees. This result may be compared 
with reference [5] where it is stated that recognition takes 
place if the positioning cf a Vander Lught filter is 
accurate to +3 degrees. In the present experiments a 
recognition spot was observed when the angular position of 
the filter was inaccurate by +43 degrees but its power was 


only 65-75% of maximua. 


Finally, it should be mentioned that throughout 
the experiments an improvement was observed due to the use 
of a liquid gate for the objects. The objects were used 
many times fcr different experiments and after each use they 
were washed in detergent to remove any oil left from the 
uSe. In the first experiments it was found that the power 
in the recognition spot decreased by about 20% when a liquid 
gate was not used. It was found, however, that this 
percentage slowly increased as the experiments progressed 
and reached 40-45% at the end of the experiments. The 
reason for this phenomencn is very likely that small 
scratches and imperpecticns develop as the objects are used 
and washed, reducing the efficiency of the object. This 
reduction in efficiency could be quite serious ina 
practical application if liquid gates are not used since a 


20% reduction may fossibly be tolerable but not 40-45%. 


6.5 Discussicn 


ee ee ee ee oe oe 


The first set cf experiments demonstrate 


convincingly the ability of the Lohman filter to identify 
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fingerprints. For each of the three filters a visual 
inspection left no doubt that the print which the filter was 
matched to was the only print which caused a recognition 
spot to appear in plane F3. The power measurements 
substantiate the above statement but it should be pointed 
out that the power measurements do not really tell the full 
Story. When matching did not occur the light distribution 
in plane P3 was fairly slowly varying whereas a sharp point 


of light was cbserved when matching occurred. 


The results of this experiment may be compared 
with the results mentioned in the introduction where 64 
Vander Lught filters were tested with a file of 10,000 
prints. Since in 65% of the cases the Vander Lught filter 
identified the correct print alone out of 10,000, it can be 
seen that while the present experiments confirm the ability 
of the Lohman filter to identify fingerprints,no specific 
conclusions as to its capability for prescreening a larger 
number of fingerprints can be reached from the small number 


of measurements. 


The only way to settle the above problem is to 
develcp criteria in terms of power in plane P3 and the 
appearance of a recognition spot. These criteria would 


have to be based on a larger number of experiments. 


The positioning of the filter in its plane is one 
of the aspects of optical filtering which frequently causes 


concern. The experimental results typified by Figure 17 
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show that the positioning accuracy required for the present 


experiments is easily achieved. 


The experimental system described can be used to 


process more than one print at the same time. If an array 


of prints, each prepared in the way previously described is 


presented as an object and one of these prints matches the 


filter then the location of the recognition spot tells which 


of the prints matches the filter (the location of the 
recognition spot moves as the correct print is moved in 


plane P1). 


As the number cf prints processed in parallel 
increases, the accuracy required in positioning increases. 
The reason fcr this is that while each print in the input 
has the same bandwidth the spatial extent increases. The 
"“space-bandwidth product" (SB product) is a measure of the 
total amount of informaticn contained in the object. All 
the information in the okject is also contained in its 


Fourier transform. The overall size of the Fourier 


transform of the object does not increase as the SB product 


increases. So, to accommodate the increased information, 
the structure of the Fourier transform becomes finer and 


thus positioning becomes mare critical. 


It is fcund experimentally that with the fresent 
system the filter can be positioned seven times more 
accurately than is required for one print. Thus if only 


the cross section of the laser beam is expanded the systen 
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is capable of processing 49 (7x7) prints in parallel. 


The experiments with angular positioning show that 
in order to maintain 90% of the maximum power in the 
recognition spot the angular positioning accuracy must be 


approximately +1.5 degrees. 


When a fingerprint is found at the scene of a 
crime the police are generally able to pinpoint its 
orientation tc +15 degrees [25]. Thus for the purposes of 
criminal investigation a print would have to be tried at 10 
different angular fositicns in order to process it 
completely. This could be done by constructing 10 
different filters and processing the print sequentially by 
these 10 filters in separate systems or separate branches of 
the same machine. Alternatively the parallel processing 
capability of the filtering system could be used to this 


end. 


The aspect of the experimental results which are 
likely to cause the most difficulty in practical use is the 
dependance on a liquid gate to insert the object into in 
order to achieve consistent and accurate measurements. To 
mount each object of a large collection in a liquid gate is 
unpractical and to mount each object just prior to 


measurement is time consuming. 


When the objects were new and unworn the power in 


the recogniticn spot was found to decrease by about 20% when 
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a liquid gate was not used for the object. However, a 
recognition spot was still observed. Thus it may ke 
possible to do without a ligiud gate and obtain useful 


results. 


It appears that in order to increase the useful 
life cf an object other cleaning methods or alternatively, 
better recording materials than the photographic film have 


to be develored. 


90 


antRnh see cakes 
owed £27 SbldeTpOIOds eds nes? 


A 


rd 


CHAPTER 7 


SUMMARY AND CONCLUSIONS 


In recent years optical spatial filtering and in 
particular optical matched filtering has found many 
practical applications and the interest in this field has 


been considerable. 


A convenient way of producing filters capable of 
controlling the spatial distribution of both amplitude and 
phase of light transmitted through them was for a long time 
the main obstacle in the field of matched filtering. This 
problem has keen greatly reduced in recent years by the 
development of methods tc produce a special type of 
holograms or masks, consisting of absorption patterns only 


Which can be used as matched filters. 


The first of these methods is an optical 
(interferometric) method due to A. Vander Lught. Many 
experiments, some of which are mentioned in the 
introduction, have been carried out using filters produced 
by this methced and the pctential and limitations of these 


filters are well documented. 


A second method of producing matched filters is 
due to A. W. Lohman. This method employs a digital 


computer and a computer guided plotter to produce a filter 
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from a mathematical knowledge of the filtering operation. 


In many applications a Vander Lught filter can be 
very easily produced. For this reason less attention has 
been paid to matched filters produced by Lohman's method 
than to the Vander Lught filters and no study of the 
performance of the Lohman filters as matched filters has 


been published. 


While the methcd of production is quite different 
for the two filters, they appear to have many properties in 
common. However, no study of these similarities of the two 


filters has been published. 


In this thesis a theoretical study of the 
transmission properties cf the two filters has been carried 
out and it has been Shown that the transmission functions of 
the two filters are indeed very similar. The similarity is 
so strong that the Lohman filter can be regarded as a 


computed binary simulaticn of the Vander Lught filter. 


The above theoretical studies indicated that in 
most filtering problems either a Vander Lught or a Lohman 
filter could be used. The choice as to which filter is 
used would depend mainly on the ease and economy of 
production. There are, however, many instances where the 
use of a Lohman filter is advantageous or a Lohman filter 


can more easily be produced than a Vander Lught filter. 


In order to study experimentally some of the 
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properties of Lohman filters, fingerprints have been used as 
test patterns in the present studies. Initially a survey 
was carried out to determine the sampling density required 
from a fingerprint to compute a Lohman filter. It was 
shown that a sampling rate of 5.5 samples/mm was sufficient 
and thus a compiete fingerprint can be represented by 
256x256 samples. In this connection the cost of producing 
a Lohman filter has been outlined and it was found that the 
cost of producing a typical filter consisting of 12&x128 
elements is approximately $25-30. Using the sampling rate 
previously determined, some filters of fingerprints were 
made according to the Lohman method. These filters were 
experimentally tested to study their properties and the 
practical difficulties in their application to fingerprint 
identification. The first experiment performed with these 
filters was designed to test the ability of the filters to 
identify fingerprints. In the experiment each of three 
filters was used to identify one print from a collection of 
eight. The results of the experiment were very encouraging 
Since in all cases a sharp recognition spot was observed 
when the print which the filter was designed to detect was 
being considered. No other print caused a visual 
recognition spot and the maximum power measured where the 
recognition spot should have appeared was only 14% of the 
power measured when recognition took place. EX periments 
were next performed to investigate the requirements of the 


experimental system. It was found that in order to obtain 
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good results in the recognition experiment the filter had to 
be accurately positioned in the appropriate plane of the 
recognition system. The accuracy required in the 
positioning normal to the optical axis of the system was 
found to be 475 microns while it was found that the filter 
and object had to be adjusted to within +1.5 degrees in the 


relative angular position. 


The experimental system used could be adjusted to 
+10 microns normal to the optical axis and an angular 
position of +10 minutes. Thus the accuracy required is 
easily achieved. It should in fact be possible to design a 
device to automatically change the filter and position it 


with the required accuracy. 


It was shown that the experimental system is 
capable of processing many prints in parailel and with the 
achievable accuracy in the present system up to 49 frints 


can be processed at the same time. 


The experimental result which was found to 
potentially cause the most difficulty was that during the 
experiments the filter and preferably also the object had to 


be mounted in liquid gates. 


It was found that if a liquid gate was not used, 
the efficiency of the reccgnition process decreased 20-45% 


depending on the conditicn of the object. 


In many applications, for example, in fingerprint 
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identificaticn it may be possible to tolerate a 20% drop in 
efficiency if other checks, such as human observation, are 
increased. However, a 45% drop in efficiency would in most 


applications be excessive. 


Since the mounting of a large collection cf 
fingerprints in liquid gates would be exceedingly cumbersome 
and expensive the use of liquid gates similar to the ones 
used in the present experiments would be disadvantageous. 

It would therefore be important to develop recording 
materials which de net reguire the use of liguid gates. 
Alternatively, it might ke possible to design more 
convenient liguid gates and develop methods for preserving 


the object in a gcod ccndition despite repeated use. 


Finally it may be pointed out that while the 
requirement of liquid gates might be a serious limitation 
due tc the large number cf objects that would likely have to 
be tested, there are many applications for Lohman filters 
where the speed with which the object (and/or filter) can be 


changed is net critical. 
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APPENDIX A 


It is required to evaluate the inverse Fourier 
transform of equation (3.11) which is reproduced here for 
reference. 
u-(n+P_)d v-nd 

m 


Rect Rect 
cd Wed 


b(x,y)=F-? Exp (j27x,u) 


| 
Brg 


(34041) 


The inverse Fourier transform of equation (3.11) 
can be obtained by Fourier transforming equaticn (3.11) and 
then changing the signs cf the coordinates in the Fourier 
plane. In Figure 5 the coordinates in plane P3 are negated 
with respect to the coordinates in plane P2 so that a direct 
Fourier transform expressed in the reflected coordinates is 
in fact the inverse Fourier transform (c.f. section 2.3). 

It is seen that the expression under the summation signs in 
equation (3.11) can be written as a product of two 
functions, one depending on u only and the other depending 
on v only. Thus the evaluation of the Fourier transform of 
equation (3.11) is reduced to evaluating the following 
transforms: 
u=(ntP) a 
r(x) =F Rect ——_—_—_—_—— |Exp {j27x,u} (A. 1) 
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t(y) =F ali (A.2) 


In order to evaluate the above equations the shift 
and similarity theorems of Fourier transforms are used. 
These theorems are stated below in terms of the coordinates 
in equations (A.1) and (A.2) but a proof can be found in 


reference [{13]. 


Shift Theorem: 
If G(u) is a ccmplex function such that 


F{G(u)}=9 (x) and a is a ccmplex constant then 
F{G (u-a) }=g (x) Exp {-j2Txa} (A.3) 


Similarity Theorem: 
If G(u) is as defined above and b is a complex 


constant then 
F{G (u/b)}=( bi xg (bx) (A. 4) 


Finally, the Fcurier transform of a Rect-function 


is required. This is given by [22] 
F {Rect (u) }=Sinc (x) (A.4) 


Combining equations (A.5) and (A.4) and changing 


variables yields 
F {Rect (u/b) }=|b]| xSinc (bx) =| b|x Sinc(t) (A.6) 


where t=bx. Now the shift theorem and (A.6) can be 


combined which with resubstitution for t results in 
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F {Rect (u/b~-a/b) =|b]| xSinc(t) Exp {-j27ta/b} 


=| b{* Sinc(bx) Exp {-j27ax} (A.7) 


This last equation can be used to solve (A.1) and (A.2). 
Consider first (A.1). In equation (A.7) let a= (nt+P a) and 
let b=cd. With these values of a and b equation (A.7) 
gives the Fourier transfcrm of the Rect-function in (A.1). 
Applying the shift theorem once more to allow for the 


exponential factor in (A.1) results in 
r(x) =cdSinc{cd (x+x..) JExp{-j2md(n+P__) (x+x )} (A. 8) 


Next in order to solve (A.2) let a=md and b=W od in eguation 


(Au) This gives the direct result 
t (y) =dW,,Sinc (yW 4) (Ae?) 


Now in order to find the Fourier transform of the sum of 
products of Rect-functions in equation (3.11) the 
contributicn due to the various Rect-functions is summed. 


This results in equation (3.12) 


b (x,y) =cd?#Sinc[ (x+x_) cd] 


p 2W ym Sine (yw, 4) Exp {-j2md[ (n+P__) (x+x_) tny }} (3.12) 
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APPENDIX B 


OMEUTER PROGRAMS 

The computational task involved in the production 
of Lohman filters is conveniently split in two. First the 
computation of the parameters which describe the filter and 


secondly the plotting of the filter. 


Programs which perform these tasks are listed on 
the following pages. Comment statements have been 
extensively used so that the programs may be easily 
understood in context with the discussions in the main body 
of the thesis. Subroutines used in the main programs are 
listed immediately following the main program, excert the 
System subroutines READ and WRITE and the Fourier 
transforming subroutine HARM. This last subroutine is 
available in the IBM Scientific Subroutine Package which is 
stored in the public file *SSPLIB in the computing system at 
Computing Services, University of Alberta. Subroutine PLOT 
is in the IBM Calcomp Sukroutine Package stored in the 


public file *ELOTLIB. 


OH 


ao Betehl eis gyeed saat w20tieq dotdy aaszyoxt er 1 
eee 2anoustare snares er td 


yood sits / dit af snobadtonelS edt do4v axadada at Beodasebau 
+78 naesyOm ara sit at ben os ttwoanwe ~akeods edt 20 
aay Hyenns ,#6spoTy atne act palvolsoe thesekboaez beteall 
teixyol odt fos ATEN bas GAAG seatiuordue weteya 

et snigdo1dhe test afd? ~~ mean ‘sntidordua pabesetannss 

et A9L4W BHD SnIuordee 3iheastve WAT add BE @ldal fads 
i5 G6far2 walfidqmo> eds at GLiee* sif¥ SEbdog ad? ‘te bezosa 
TOG sufsvesrdta sh¢460Lh Fo qitatavind jsesivese pais uqaoo 
ad3 at feielte ape#989 9als¥orlue qaootnd mar ad? ai et 


CVO Qa @Q Oi@r@ (Ve) (oe) er@) Oe (alternate! 


Aaa Q OO 


77 


10 


102 


THE FOLLOWING IS A FORTRAN SOURCE PROGRAM TO FOURIER 
TRANSFORM AN ARRAY OF NUMERICAL PICTURE SAMPLES AND 
COMPUTE THE PARAMETERS WHICH DESCRIBE A LOHMAN FILTER 


DIMENSION A(256, 128) »M(3) , INV (4096) ,S(4096) , VA (32768) 
INTEGER*2 VP(16384) ,LEN 

COMPLEX CA(128,128),C (64,64) , CMPLX 

EQUIVALENCE (VA,CA), (VA,A) 


READ INPUT DATA FRCM A FILE 


DO 77 I=1,128 
K= (I-1) *128+1 
CALL READ (VP (K) ,LEN,0, LNR, 1) 


THE INEUT FILE CONTAINS INTEGERS 0 TO 63, CONVERT 
THESE TO REALS 0.0 TO 1.0 AND EQATE TO THE REAL PART 
OF THE COMPLEX NUMEEFRS TO BE INPUTTED TO HARM. 


DO 10 I=1, 16384 
VA (2*I-1) =FLOAT (VP (I)) /63. 


SET INPUT PARAMETERS FOR SUBROUTINES HARM AND SFLE 


N1=128 
N2=128 
M1=64 
M2=64 
M(1)=7 
M(2)=7 
M (3) =0 


REARRANGE INPUT FILE TO A SUITABLE FORN FOR HARM, CALL 
HARM TC PERFORM FOURIER TRANSFORMATION AND REARRANGE 
OUTPUT FILE 


CALL SFLE(CA,N1,N2,M1,M2,C) 
CALL HARM(A,M,INV,S,1, IFERR) 
CALL SFLE(CA,N1,N2,M1,M2,C) 


CONVERT THE REAL AND IMAGINARY PARTS OUTPUTTELD FROM 
HARM TC AMPLITUDE AND PHASE NOTATION MAKING SURE THAT 
IF BOTH REAL AND IMAGINARY PARTS ARE ZERC, CONVERSION 
IS NOT ATTEMPTED AND AMPLITUDE SET TO ZERO. 


PI =3. 14159 

DO 21/1=1)128 

DO 21 J=1,128 
RP=REAL (CA (I,J) ) 
CP=AIMAG (CA (I,J)) 
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IF (RP.EQ.0..AND.CP.EQ.0.)GO TO 21 
AMP=CAES (CA(I,J) ) 
PHA SE=ATAN2 (CP,RP) / (2. *P1) 


THE ABOVE STATEMENT NORMALIZES THE PHASE TO VALUES 
BETWEEN -1/2 AND 1/2. K IN EQN. (3.23) HAS BEEN 
TAKEN AS 1.0. 

NEXT STORE AMPLITUDE AND PHASE IN CA(I,J) AND GO 
BACK TO CONVERT NEXT NUMBER. 


CA (I,J) =CMPLX (AMP, FHASE) 
CONTINUE 


NEXT NCRMALISE AMPLITUDES TO VALUES LT.OR.EQ. 1.0, BY 
FINDING THE LARGEST AND DEVIDING ALL BY THAT VALUE. 


BIGAMP=0. 

DO 2291517128 

DO »22*3=1,.128 

IF (REAL (CA (I,J) «GT. BIGAMP) BIGAMP=REAL (CA (I,J) ) 

DO 23 I=1,128 

DO 23 J=1,128 

CA (I,J) =CMPLX (REAL (CA (I,J) /BIGAMP, AIMAG (CA (I, J))) 


FINALLY STORE THE FARAMETERS IN A FILE WHERE THEY CAN 
BE USEC BY THE PLOTTING PROGRAM. 


LEN=40 96 

do 20) j=1,32 

K= (J-1) *1024+1 

CALL WRITE (VA(K) ,LEN,0,LNR,2) 
STOP 

END 


SUBROUTINE SFLE INTERCHANGES 1ST AND 3RD QUADRANT AND 
2ND AND 4TH QUADRANTS OF THE INPUT/OUTPUT MATRIX OF 
SUBROUTINE HARM 


SUBROUTINE SFLE(CA,N1,N2,M1,M2,C) 
COMPLEX CA(N1,N2) ,C(M1,M2) 


STORE ZND QUADRANT IN C(1I,J) 

DO 1 I=1,M1 

DO 1 J=1,M2 

C (I,J) =CA (I,J) 

MOVE 4TH QUADRANT INTO 2ND QUADRANT 
DO 2 I=1,4M1 

DO 2 J=1,M2 

CA (I,J) =CA (I+M1,J+ M2) 


MOVE C(I,J) INTO 47TH QUADRANT 
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DO 3 I=1,M1 
DO 3 J=1;N2 
CA (I+M1,J+M2) =C (I,J) 


STORE 1ST QUADRANT IN C(I,J) 


DO 4 I=1,M1 
DO 4 J=1,M2 
C (I,J) =CA (I, J+M2) 


MOVE 3RD QUADRANT INTO 1ST QUADRANT 


DO 5 I=1,M1 

DO 5 J=1,M2 

CA (I,J+M2) =CA (I+M1,J) 
DO 6 I=1,M1 

DO 6 J=1,M2 

CA (I+M1,J) =C (I,J) 
RETURN 

END 
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THE FOLLOWING IS A FROGRAM TO PLOT THE ENLARGED 
ORIGINAL OF A LOHMAN FILTER. THE PROGRAM MAKES USE OF 
ONE SUBROUTINE, PLCL 


DIMENSION VA (32768) , WORK (2048) 
COMPLEX CA(128, 128) 

INTEGER*2 LEN 

EQUIVALENCE (VA,CA) 


SET THE SIZE OF EACH CELL (DIMENSION OF CELL= U X UD) 
AND THE THICKNESS CF PEN TO BE USED FOR PLOTTING 


U=.22 
PEN=.0 14 


READ FILTER PARAMETERS FROM FILE WRITTEN BY 
PROGRAM PREVIOUSLY LISTED AND SCALE THE AMPLITUDES 
IF REQUIRED 

RMULT=4. 

CHOP=2. 

DO 88 I=1,32767,2 

VA (I)=VA(I) *RMULT 

IF (VA(1) .GT.CHOP) VA(I) =CHOP 


INITIALISE PLOTROUTINES AT A LOW ADDRESS IN CCRE AND 
SET ORIGIN AT TOP LEFT HAND CCRNER OF FILTER, ALLOWING 
FOR A .5 INCH MARGIN. 


CALL PLOTS (WORK,8192) 
CALL PLOT(.5,128.*U,-3) 


INITIALISE A COUNTER FOR THE COLUMNS PLOTTED AND SX, TH 
DISPLACEMENT FROM THE ORIGIN IN THE X-DIRECTICN OF THE 
LEFT TCP CORNER OF THE CELL TO BE PLOTTED. THEN PLOT 
THE FIRST COLUMN. SUBROUTINE PLCL PLOTS THE APERTURE 
WHEN GIVEN (SX,SY), THE COORDINATE OF THE TOP LEFT 
CORNER OF THE CELL AND THE AMPLITUDE AND PHASE 
PARAMETERS 


J=J+1 

SX=SX+U 

DO910, 1=4,128 

SY=~FLOAT (I) *U 

AMP=REAL(CA(I,J)) 
PHASE=AIMAG (CA (I,J)) 

CALL PLCL(SX,SY,AME, PHASE, U, PEN) 
CONTINUE 


PLOT THE NEXT COLUMN FROM BOTTOM UPWARDS 
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J=J+1 

SX=SX+U 

DO 20 I=1,128 

SY=-FLOAT (129-1) *U 

AMP=REAL (CA(129-I,J)) 
PHASE=AIMAG (CA(129-I,J)) 

CALL PLCL(SX,SY,AME,PHASE,U,PEN) 
CONTINUE 


IF ALL COLUMNS HAVE BEEN PLOTTED STOP, OTHERWISE 
CONTINUE WITH NEXT COLUMN. 


IF(J.L1.128)GO TO 9 
STOP 
END 


SOUBROUTINE PLCL 


SUBROUTINE PLCL DECIDES WHICH TYPE OF APERTURE IS 
REQUIRED FOR EACH CELL AND PLOTS IT. 


SOBROUTINE PLCL (SX,SY, AMP, PHASE, U, PEN) 
IF AREA REQUIRED IS LESS THAN 1 PENSTROKE GO TO 10 
IF (AMP.LT.2.*PEN/U)GO TO 10 


EVALUATE X1,X2,Y1,Y2 WHICH DETERMINE CORNERS OF 
APERTURE 


Y1=SY+PEN 

Y2=SY+U 

DX=.5*U*AMP 
X1=SX+ FHAS E*U-. 5*DX+PEN 
X2=SX+ PHASE*U+.5*DX 


DRAW THE OUTLINE FRAME OF THE APERTURE 


CALL PLOT(X1,Y1,3) 
CALL PLCT (X1,Y2, 2) 
CALL PLOT (X2,Y2,2) 
CALL PLOT (X2,Y1, 2) 
CALL PLOT(X1,¥1,2) 


IF WIDTH OF APERTURE IS LT.OR.EQ. TO 2 PENSTRCKES 

THE PLOTTING IS COMPLETED. OTHERWISE FILL THE FRAME IN 
WITH VERTICAL LINES SPACED 1 PENWIDTH APART, CHECKING 
AFTER EACH LINE IF PLOTTING IS COMPLELED. 


X1=X1+ EEN 

IF (X1.G2.X2) X1=X2 

CALL PLOT (X1,Y1, 2) 
CALL PLOT (K1,Y2,2) 

IF (X1. £Q.X2)GO TO 1000 


q2Lwaauro ,q0Te vstous sai ty 


Be | 


2I s807TaH i) oa a5 aw suk ici 
alee 1 mosd with 


i Riieieehpclehay: inp i 


OF oF de IHOSRa4SS F WAHT REEL CT castagen Anam Ae 
or or 98 NESEY Sean ae ; 


= 
10 2AaWHOD FHLHGETNO HOTUW Ex Sr .tt tebe hiaaas 


sanoRTenad. sD .Om, 
41 (SARt. TAT 1115 S2tua 


QiQr@r@s@iore 


10 


1000 


107 


X1=X1+PEN 

IF (X1. GT.X2) X1=X2 
CALL PLOT (X1,Y2,2) 
CALISPEOT(L 1, 0172) 
IF(X1.EQ.X2)GO TO 1000 
GO TO 1 


THE FOLLOWING STATEMENTS PLOT APERTURES WITH AREA 

LESS THAN 1 PENSTRCKE. FIRST CHECK IF AREA 

REQUIRED IS LESS THAN THAT OF THE PENTIP, IF SO 
NOTHING IS PLOTTED. OTHERWISE DRAW A LINE OF THE LENGT 
REQUIRED TO PLOT TEE DESIRED AREA. 


PI=3.14159 

CH=. 5*AMP* (U**2) - ( (PEN/2. ) **2) ¥PI 
IF(CH.LT.0.)GO TO 1000 
DY=CH/PEN 
X1=SX+ PHAS EXU4+PEN 

Y 1=SY+PEN+ (U-DY) /2. 
Y2=SY+FEN+ (U+DY) /2. 
CALL PLOT (X1,Y1, 3) 
CALL PLOT(X1,Y2,2) 
RETURN 

END 
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